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HAUPTAUFSÄTZE 
Allgemeine Schalentheorie I 


Von H. Neuber, Dresden | 


‚Ausgehend von den Grundgleichungen des elastischen Kontinuums wird mittels des absoluten Dif- 

. ferentialkalküls eine allgemeine Theorie für Schalen beliebiger Form und Berandung unter Berücksichtigung 

der Veränderlichkeit der Wandstärke abgeleitet. Durch Binführung konvarianter Raumableitungen. des 

Spannungstensors an Stelle der sonst üblichen Schalenlängskräfte und -momente wird eine beliebig hohe 
Approximation an den dreidimensionalen Zustand ermöglicht. 


Based on the Ihree dimensional equations of elastic media a general theory of shells of any shape and 
boundary with regard to variable thickness is developed by means of the absolute differential caleulus. Using 
covariant three-dimensional derivations of sireess tensor instead of the usual forces and moments the theory 
allows an approximation of any degree to the threedimensional state of stress. 


RE Au moyen du calcul differentiel absolu une theorie generale des Eouelles d’une forme quelconque et 
d’epaisseur variable est developpee ü la base des Equations trois-dimensionelles elastiques. Par introduction 
des composuntes derivees covariantes du tenseur trois-dimensionel au lieu des forces et moments ordinaires 


la theorie permel une approximation de queleongue degree a V’&lat Elastique exact. 


Mexona u3 OCHOBHLX ypasHenmä TeOpum YHPYIOCTU CTPoHuTcCA B IIPOHSBONBHEIX 

KPHBONNHEÄHBIX KOOPAHHATAX O0IMAA Teopus ODONOYeR I106OU POPMELH IIEPeMEeHHON TONINHHEL. 

: Ilpn momomm. koBapnautHoä TpexMmepHoi IPOU3BOnHOH TeHsopa HANPA:KeHHÄ, KOTopan 

Deperca BMEeCTO CHA H MOMEHTOB B OOBIYHOH TEOPHH, MOSKer ÖBITB NOCTUTHYTO IPnÖ,n:keHne 
Br 10004 .CTeEIeHH TOYHOCTH K TPeXMePHOMY HANPSGKEHHOMY COCTOAHNIO. 


1. Einführung N 

Die Entwicklung der Technik hat auf den verschiedensten Gebieten zu neuartigen Bau- 
teilen mannigfacher Form geführt, welche hinsichtlich ihrer Dimensionierung an den Ingenieur 
stets wachsende Anforderungen stellen. Dies gilt insbesondere für Bauteile mit verhältnismäßig 

geringer Wandstärke, die sog. „Schalen“. Die klassische Schalentheorie bezog sich vornehmlich 
auf den Behälterbau und lieferte fast ausschließlich Berechnungsgrundlagen für Rotations- 
schalen!)2)®). Eine später von Trefftz*) gegebene Ableitung der Schalenbiegungs-Glei- 
chungen auf Grund des Castiglianoschen Prinzips beschränkte sich auf die Krümmungslinien 
der Schalenmittelfläche. 

Bei Behandlung von Schalen beliebiger Form und Berandung macht sich die Beschränkung 
der Theorie auf spezielle Koordinantensysteme nachteilig bemerkbar; denn — wie bei den 
meisten technisch wichtigen Problemen der Elastizitätstheorie —ist auch bei der Schale die 
Einführung eines günstigen Koordinatensystems für die erfolgreiche Durchführung der Rech- 
nung von großem Vorteil. So entstand das Bedürfnis nach Erweiterung der Schalentheorie auf 
beliebige Koordinaten. Angeregt durch neuzeitliche Festigkeitsprobleme des Flugzeug- und 
Motorenbaues, entwickelte der Verfasser im Jahre 1939 ein dieser Forderung angepaßtes Ver- 
fahren für Schalen. allgemeiner Form und Berandung 5). Zufolge der Nicht-Orthogonalität 
eines beliebigen Koordinatensystems außerhalb der Schalenmittelfläche ergab sich dabei die 
Notwendigkeit der Verwendung schiefwinkliger Koordinaten. Durch weitgehende Anwendung 
der Tensorrechnung wurde diese Theorie später in eine dem absoluten Differentialkalkül an- 
gepaßte Form gebracht, welche den Gegenstand des vorliegenden Aufsatzes bildet. Ihr Kern- 
punkt ist nicht mehr — wie in der klassischen Schalentheorie — die Bestimmung der Schalen- 
schnittkräfte und -momente. Vielmehr wird die Herleitung der Grundgleichungen aus der 
allgemeinen Kontinuumsmechanik durch Integration der nach Potenzen des Abstandes von 
der Schalenmittelfläche entwickelten Divergenz des Spannungstensors gewonnen, wobei nur 


1) Love, Lehrbuch der Elastizität, deutsch von A. Timpe, Leipzig und Berlin 1907, $ 326. 

2) Meißner, Phys. Zeitschr. 14 (1913), 8.343 ff.; Vjschr. naturforsch. Ges. Zürich 60 (1915), 8.23, 
3) Flügge, Statik u. Dynamik der-Schalen, Berlin 1935. 

4, Trefftz, Z. angew. Math. Mech. 15 (1935), S. 101#f. Wr 

5) Neuber, Schiefwinklige Schalentheorie, Ber. d. Fachtagung d. Lilienthal-Ges., Aachen 1939. 


7 


Re 


rd 


ae , eh ” \ 2 5 


solche Größen als Flächentensoren benutzt werden, welche onen! rer 
in der Kontinuumsbetrachtung auftreten. Der so ; Rechnungsweg tet eine 
beliebig hohe KAenany an den a a E pannungszustand und führt zu 
re leichungen, deren Aufbau davon unabhängig ist, ob der verwendete Spannungstensor 
ara: gemischt oder kovariant ist. Hierin liegt der grundsätzliche Unterschied der 
stellung gegenüber dem Ansatz von Lo vet), sowie den inzwischen erschienenen, ebe 
auf dem Differentialkalkül fußenden Arbeiten von Reutter‘)und Deuker’?), vo 3 
chen vorliegende Th&orie auch im Endergebnis abweicht. Ferner geht die hier gegebene Dar- 
stellung auch insofern über die bisherige Literatur des Gegenstandes hinaus, als die Verä Da 
keit der Wandstärke exakt Berücksichtigung findet und sowohl die Oberflächen-, wie "die 
Randbedingungen in allgemeiner Form angegeben werden. NN 
Aus drucktechnischen Gründen erscheint die Arbeit in zwei Teilen. Der erste bringt die 
geometrischen, Grundlagen, die Entwicklung des Spannungstensors und seiner Divergenz, “ 
sowie die Gleichgewichtsbedingungen. Der im nächsten Heft dieser Zeitschrift erscheinende 
zweite Teil enthält die Entwicklung der Formänderungsbeziehungen und der allgemeinen Rand- 
bedingungen. A; | 
; 2. Raumgeometrische Grundlagen 
Es mögen hier dieselben Bezeichnungen Verwendung finden, welche bereits in der Arbeit - 
des Verfassers über die Grundgleichungen der elastischen Stabilität ®) mit Vorteil benutzt 
worden sind. Danach sind <= 2, y=2, z=&%, allgemein %*, x, zP oder 22 kartesische 
Koordinaten, während w= «4, v= a2, w= «°, allgemein x*, =%, =”, x oder =°, die krumm- 
linigen Koordinaten kennzeichnen. Die kartesischen Koordinaten können als differenzierbare 
Funktionen der.neuen Koordinaten «, v, w vorausgesetzt werden. Für die Ableitung gilt die 
Abkürzung 


oxk 
Ba = ck wi Weiner rg ER DE a (1). 


Jedem Vektor mit dem Komponenten A, — 4% in kartesischen Koordinaten sind in den krumm- 
linigen Koordinaten die Komponenten ; 


3 | 
A= 2 4: = 4A: a Fe PR a (2) 


er; - f 
ht nn we un un ee ne = 


ni zugeordnet, wobei das Summationszeichen weggelassen werden kann, wenn für die weitere 
Rechnung vereinbart wird, daß stets über solche Indizes zu summieren ist, welche zugleich 
oben und unten auftreten. Für die Transformation eines Tensors zweiter Stufe gilt 


s Bau — Br & el rue ken 9 oe ie ee (3) 
Der metrische Tensor wird 
3 
RP HE RE ER EN 
Au 2 Pa a a Eee (4) 


mit den kartesischschen Komponenten 


mtl, für k=1, 


We ı 8). 


Hierbei ist öf das Kroneckersche Symbol. Weiter haben die Größen c;, welche als Unter- 


determinanten der Größen c*, dividiert durch |e%| gewonnen werden, ebenfalls Tensorcharakter 
und es gilt 


8 
k oA —. Ak kei — A 
sol, da— Ök, & EI RER (6). 


Wird der Wert der Determinante des metrischen Tensors mit gan 
den Wert der Determinante der Transformationsgrößen 


al='77, Rune (7) 


Ferner wird der schiefsymmetrische Tensor 


| = 9 bezeichnet, so folgt für 


1 für. Rd, m = 1,2, 8; 2a US Ne 
Elm = Ein = !— 1 für ki, m=1,.8, ar ee I> NS: 5 ae (8 
0 für k=1 oder Z=m oder k= m HARE 


‘) Reutter, Z. angew. Math. Mech. 22 (1942), 8. 87 ff, 
‘) Deuker, Z. angew. Math. Mech 23 (1943), 8. S81f£. 
) Neuber, Z. angew. Math. Mech. 23 (1943), S. 321 ff. 


göinischten en eines Tensors, dessen kovariante Komponenten durch Ber Kr; 
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als Bene Komponenten once END j ae 
Nie ei: einer ne) 2) entsteht der Vektor. 
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Hierbei und i im folgenden Er die Differentation durch einen Rereslen Deppkiskrieh, De re 


een wobei nach der Koordinate zu differenzieren ist, deren Index hinter dem Strich unten _ A 


. angegeben ist. Bei Differentiation eines Vektors A; in kartesischen Koordinaten a der E 
a Tensor Mi Erde & ® een Na 

> \ { 9 ° a L 7 . a ar “ Y7 7 
a ee = Ale are er m. RB NEE R u. 


mit de krummiinigen Komponenten, | ; WIR, PR“. 


Be 


Hierbei ist 


Ve I „ac . gve oda we! ON u 
De zn km DR, Er BER a ER 5 S 
eine Abkürzung, welche dem Christoffel schen Dreizeiger- Symbol zweiter Art entspricht. 
Der vertikale Doppelstrich kennzeichnet daher bei Differentiation von Tensoren den. voll- 
ständigen Ausdruck der Tensorkomponenten in allgemeinen Koordinaten, d.h. einschließlich 
der mit den Christoffel-Symbolen behafteten Zusatzglieder. Die so entstehende, als 
kovariante Raumdifferentiation oder hier kurz „Raumableitung‘“ bezeichnete Operation 
liefert daher stets neue Tensoren. Bei Anwendung auf einen Tensor zweiter Stufe ergeben sich 
u Rechenregeln: 


2 B Ä 
x Bul=75 en: „— Tv Bio» 
Fa, 
Bl = u + Ivy. Bu—T% 2 EN a LER 
2 
Br,=7 ve BB" 


Für die zugehörigen N te wird im folgenden die abgekürzte Schreib- 
weise 


Bu=g@eBule, Buf=weBil, 3#=WeB"l,..... (N 


verwendet. 

Bei Anwendung auf den metrischen Tensor ergibt sich infolge des Verschwindens der 
zugehörigen Komponenten in kartesischen Koordinaten, daß seine kovarianten Ableitungen 
in jedem Koordinaten-System verschwinden, d.h. es gelten die Ricci- Gleichungen °): 


OR 


9) Vgl. hierzu, sowie bezüglich weiterer Einzelheiten: Levi-Civita, Der absolute Ditferential- 
kalkül, Grundl. d. math. Wiss., Bd. 28, Berlin 1928. 
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BA; Fa Binap Bi he durch w= 
RR inition. Kon w a, senkrechter Abstand von ihr ge 


Pt EN Er U 


2 (Die HAGER, Koofäklsien "der, Schalenmittelfläche, ee n 
genannt wird, sind hierbei ausschließlich Funktionen von 7 


mit N% bezeichneten kartesischen Komponenten des N ktors der 

die weitere Rechnung ist deshalb. die Einführung unterschiedlicher | 
dizes 1 und 2 erforderlich. Die augehörigen allgemeinen ne ac. # ” 
Dann ist: BR 


we 4. 
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Re u 
die IE RIEER von ® * Järgs der Fläche, mithin e ein Tangentenvektor der Fläche. Flierane wo 7; 
sich der Re es Tensor der Fläche zu: 


Sie gestatten ihrerseits die Berechnung der kontravarianten Komponenten des Tangeritenvektors: R ' u 


/ 02 a U un RR 
in Übereinstimmung mit der auch auf der Fläche geltenden Bedingung 
PER A AL EN Ve SR 


Das Kronecker- Symbol hat mithin, wie im Raum, so auch auf der Fläche Tensoreigen- 
schaft. 


Für die Differentiation einer Invariante / auf der Fläche gilt 14) 


0 


RA : 


Hierbei und im folgenden wird die kovariante Flächendifferentiation durch einen einfachen 
Vertikalstrich gekennzeichnet. 


Die Flächenkomponenten eines Raumvektors A, sind | 
A, nA tn u RR 
Für die Flächenableitung dieses Flächenvektors gilt 


Ga 
Alı= 3 TA Rn RER 


Hierin entsprechen die Größen 


vr 0% 70h 06, RR. aus | dans dann 
&ß da Dar) FpTE + ROT Ye 3 . +(30) 


den Gristoffel-Symbolen zweiter Art der Fläche. Für die Flächenableitung von Flächen- 


tensoren gelten entsprechende Regeln. Die Flächenableitung des metrischen Tensors der Fläche 
verschwindet identisch (Ric ci): 


Re =0, wo, er 


ird ; - u 
AR u 
gesetzt, so folgen als kontravariante Komponenten a*® des metrischen Tensors, entsprechend SE 
der RBuRE ap aPr = ö: En 
5 = amt unann ZT, 0, (28), 
BERN, a a 


ıngen wird, mögen die aften des BEN re 
R » w en ne und nmerıng) he = seine, haspepschen Ba Be 
De N ® 


BEN, 8 N er ze F ACH m&=o Se Ne ve . ER Bere 2 IR RiPR k (09 


ER . REN. l er vi { MM REN ee BE 
Es Diese  Beingungen sind erfüllt, wenn. N, aus S Mh DE ER a. 


X HR { N= md a ir CE , ER . 2 ey 7 Kr 


B bestimmt wird. Mit Einführung Se este Flächentensors. EL; 
, | Yafürao—=1,ß=2; De, i Wr ee 
Be a ne 

Ofüra=P- Ba | WE@T- N 
mit den kontravavarianten ne 


für et Dee 2: 


Pr. 

as) ot 02.5 ler: 8 Re ER ER ey 5 Bi 
% En Le er 2 ee 
kann GI. (34) in der Form BR re | 98 

| | 1 ERREN N > 
N= 5 me Kor EEE RE RN RER ER . (37) } x 


geschrieben werden, welche hinsichtlich der Summation über die auftretenden Flächenindizes 
erkennen läßt, daß der Normaleinheitsvektor eine Invariante der Fläche darstellt. ‚Ferner ist 
leicht ersichtlich, daß Gl. (32) erfüllt ist. Weiter folgt für-das Längenquadrat 


| N,N* - RER re TEEN 
Mit Hilfe der Rechenregeln | 
me OL — MO, ee eh Kerr...) 
läßt sich, hieraus leicht bestätigen, daß Gl. (33) erfüllt ist. 
Die Flächenableitung von Gl. (33) liefert | 
IN EN EEE a BE RR ARE (40). 
Der Vergleich mit Gl. (32) zeigt, daß ger Vektor N#|, ebenfalls tangential zur Fläche gerichtet 
ist und daher in linearem Zusammenhang mit dem Tangentenvektor c* stehen muß, d. h. es kann 
| I Re A ee a Br EEE 


gesetzt werden. Der hierbei auftretende Flächentensor b} trägt den Krümmungseigenschaften 
der Fläche Rechnung. Sind o, und o, die Hauptkrümmungsradien der Fläche, so ergibt AR 
unter Umgehung einer Zwischenrechnung für die mittlere Krümmung 


| RL A, a 
E; B-a(:t.)=5% RE Rue 17:42) 
$ und für das als „Gausssches Krümmungsmaß“ bekannte Produkt der Haupikrümmungen 
A der Ausdruck 

ee 0 A TE (43) 
| | EN Een 

wobei 

4 Bene NEE RE ln ee Ne ER RN (44) 
; gesetzt ist !P). 


10) Vgl. hierzu, sowie hinsichtlich weiterer Beaipeben, der Flächentheorie: B laschke, Diffe- 
rentinlgeomeitie I, Berlin 1945, 
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Bi | Is=0, 95-1 te 
Die zugehörigen kontravarianten Komponenten ergeben sich bei Beachtung der Bedingung Be 


2 7 


a & Ian 9"” = 64 in der Form: . u 
= ER AR +3W ch +..., I \ BE a Br 
>. vi Für jene Christoffe oki, bei welchen der Index 3 auftritt, ergibt sich 

1 BERL s ___ 1 29uB 

e) a | RR I ee 

Be PEN ai a 51). 
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Ts=0, Ta=0, IT; = 0 


Bei der kovarianten Differentiation der Flächenkomponenten von Raumtensoren erweist sich 
die Einführung der unterschiedlichen Bezeichnung für Raum- und Flächendifferentiation als 
besonders vorteilhaft. Bei Anwendung auf die Flächenkomponenten D, seines Raumtensors D;,, 
an der Stellew=0 folgt z.B.: 


Derllue D De Ds AT 2 Das (52). 


Bei Anwendung auf den Vektor cl} ergeben sich wichtige Gesetzmäkigkeiten der Flächen- 
theorie. Aus dem Verschwinden seiner Raumleitung folgt nämlich 


ok 
(da = a en (53) 
oder mit- Bezug auf Gl. (51) 


ie Ob a Bup NR ea (54) ' 

In gleicher Weise folgt für ck: h h 

3 (alldw=o = N’ —(Tso)o= 0%, N NEE (55), i 
‘d.h. in Übereinstimmung mit Gl. (41): 18 
k 0 . i 

Nas bb... 0 (56). 4 


Wird Gl. (54) nochmals differenziert, so ergibt sich für den Unterschied zweier Ableitungen bei 
Vertauschung der beiden letzten Indizes mit Berücksichtigung von Gl. (56): 


6a 0» — Ca |ya= (b (da 8l» — bay|) N + (by, b5 — bag b}. >) 05 ET SST. (57). 

Andererseits kann die linke Seite der Gleichung auch mit Hilfe der ‚Regeln der kovarianten 
Flächendifferentiation umgeformt werden, indem von der Identität ; 
ler 0, dw | 
rare Baer.) 


Nee sch Eurtcher die Bedehungen von 6: aus s 
N een mit Are, RE ET RE NEE 
ER a era brnla Be 7) Era a ee 
laß der sr @ (59). enthaltene Keee vierter Stufe Bean Kamm 
Ben Tensor U RER LE HTENE 
ER NEFZE #2 Rap: =“ ba ba baby u DR or en * NS Rn 
bereinstimmt. )ies ist der Ri iemann- % h Tr is t oT, e 1-Tensor der: Fläche, ‚welcher b 
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Re: "T,=1,—ubl, u? Bon. re \s ern 2. 
Banlt sind sämtliche Hiltsmitte) erörtert, welche erforderlich sind, um in der näheren .; 
der Fläche die Raumableitung der Komponenten eines Raumtensors bestimmen, bzw. auf 
nn. n Flächenableitungen zurückführen zu können. Zur Behandlung des Schalenproblems ist ferner 
die Entwicklung von Tensoren nach der Koordinate #° = w erforderlich. Bei Verwendung der 
2 .Taylor-Reihe erhält man hierbei zunächst partielle Differentialquotienten nach w. Im 
u ' Sinne der bereits in der Einführung erwähnten Zurückführung aller Spannungsgrößen auf Ko 
E: 5 Raumtensoren, welche in vorliegender Arbeit als Richtlinie zur Gewinnung eindeutiger End- 
 gleichungen konsequent eingehalten wird, ergibt sich die Aufgabe, die auftretenden Differential- DR 
Be ‚quotienten nach ® durch Raumableitungen zu ersetzen. Die hierfür erforderlichen Beziehungen 
' sind nachstehend am Beispiel eines Raumvektors 4, sowie der symmetrischen Raumtensoren 
000 Ban — BA, (%uv — QuRV und Dinve — Duive — DAxev mit Bezug auf die Fläche w = 0 wieder- 
0... gegeben, wobei ein besonderer Hinweis betreif w= 0..der Kürze halber weggelassen ist: 


R Ally — 42], ER EEEN HT DR | | Er 
All, — 42, +5,..48, | | lin! 
En: 6137 ‚4lly= = Ayo +b,. 4° lö + bs Ara Hbuledt 3, Aa 
ee ER 94° & “Ad, a e 
e ae + | RE a ER 
E.:\ FAT _ gell, 288 49, 20 40 un a 
E Ga Allen +2 Alla HA nen 
043 02.43 ! [} E 
le: alle N 


2, Ben, er b} B3# — bi Be3, 
B*®||,= B°®|, + by B*’ — b}, B®, 
B®||, = B®|, +2 b,,.BP®, | (65), 
B®||, (Er B®|ys +2 2b, Bel; +2 Dao b° ae + 2 b,ols 8 +2 Re bir DIR 
— 20,5 B®— b,; B® ||, 


ri el +53 BrP +58 Ber, 
dw 
EB ln #20; Bil Hau Blad2e But Ha m Ham, 
ee Bra], + bs Br0 (66), 
* w 
22. Be3 - = 
En — Böll, +2 55 Bräll, +20: 8°, 
Br { i 
= B%ll,;; +3 b2 Bel, +6 B’?||; 2.08 b) BY3 


au 


=. je vorn 
Or], — 01.000 +,1[0808 
en 


By "Wird al Spannungstensor der Schale mit Se bezeichnet, so il 
Schalenmitteltläche die Taylor-Entwicklung 


u = [Pr ]u=o $ w je (8% ). -0 +2 u. 


Die hierin auftretenden naeh Differentialquotienten Hack «2 kind 

(66) durch Raumableitungen zu ersetzen, wobei im folgenden der Kür 

w= 0 weggelassen und dafür der Index 0 angebracht ist. Damit die im 
zu formulierenden Gleichgewichtsbedingungen der Schale alle Glieder bis einschließlich rail sr 
ständig enthalten, muß die Entwicklung der Komponenten 8*# bis w®, der u ee 
und 8% bis w® durchgeführt werden, Mit Benützung der Gl. (66) ergeben sich dann für rücken 
‚ponenten des Spannungstensors ye Schale N ‚Ausdrücke: N 


ae +28 |, +) 


a a a 9 Hören) +} 211. 
0 
+ 20,897 rag Hass 2 


“ 
a 


5 un, 5 G, “ ki 
IR sms tee nz Haus raume] lan 
mi... ee TERN 


een 


333 


An der Ober- und Unterfläche der Schale stehen die Spannungskomponenten der Schale mit 
den äußeren Kräften im Gleichgewicht. Im folgenden Abschnitt werden die zugehörigen Be- 
a. unter Berücksichtigung der Veränderlichkeit der Wandstärke angegeben. 


5. Die Oberflächenbedingungen der allgemeinen Schale veränderlicher Wandstärke 
Für die kartesischen Koordinaten der Ober- und Unterfläche der Schale gilt mit Bezug 


auf Gl, (20), wenn die Wandstärke der Schale mit % bezeichnet, also w = -+ Z u wird: 


s wi 
Demnach folgt für den a I der Ober- bzw. Unterfläche: 


*, 0 
KV=-Kr- “ ng 5 te Im 1 RE EN (72). 


Für den Normaleinheitsveklor der Ober- bzw. Unterfläche gelten die Bedingungen (Orthogonali- 
tät und Normierung): | 


NO, N, N Te ah 


Hieraus ergibt sich unter Vernachlässigung höherer Potenzen der Wandstärke bei Umgehung 
der Zwischenrechnung: 
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N = (1— 2A) wer — Pu 2 


auf, Werden die Komhente der an den Öber- und Unterfläche Nur Schale eier Anbenkti x 
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” ee Spnmungslomonaten auf die Ober“ baw. nteräche der sch: 


sh Ka AR N “Hr Nr } 
“ hteir =; a er ee . A 3, Fey a las un A 

mit den. Komponenten RE, RE OR VICKTE ei $ AL, N. 
er 1 s ur) R Ale, ar, Ph 


De =: +5 Ei SERTUAS Gas — 1 nn np. ae BR i = Baer A (76). E 


‚Kräfte mit k“ + I bezeichnet und auf die Fläche w — konst. bezogen, auf welcher der be- ur 


22 ; Bere Obertischenpun liegt, so ergeben sich folgende Randbedingungen: ne ER 


Gar) a = dos (br ae IREZER: m) x ! | 


Nach Division durch 995 und Einsetzen der Spannungskomponenten aus Gl. (70) lem hier E 
mit Bezug auf die Gleichungen 76) im Rahmen der Vernachlässigung höherer FOR sat va 
ie und ihrer Me . 


+ be 8 +80 +4 


Er 5 
ee 2) 
A 2. 


„rast, +0 Ser 


RR Hase +63) 2] £ 


3er mpterjerbe sepjterla a 
8). 

$ = be Rn 4a, Fa. Hua BEN Be 
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33 
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Hierbei bilden die Glieder mit festem Vorzeichen und jene mit wechselndem Vorzeichen zwei 


“ getrennte Gleichungssysteme. Die Auflösung nach Dr 883 El und 8%|| führt bei wieder- 


0 0 
holter gegenseitiger Substitution im Rahmen der Vernachlässigung höherer Potenzen der Wand- 
stärke und ihrer Ableitungen auf folgende Ausdrücke: 


gr BZ | Fler Hs, 
0) 
sl, = 0% une ae 
333 
Hell rs) Has, 1...) 


AR u 3 dr 2 7% &y 
83 k En Hr za [2! En ib 
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S” ER hl rn fa Ajojge? 
I, ar Beh |, 


+38) 

333 ne 

Die links stehenden Größen sind damit auf die Komponenten der äußeren Belastungskräfte der 
Schale, sowie auf weitere Komponenten des Spannungstensors und seiner Raumableitungen 
nach w zurückgeführt. Die beiden ersten dieser Gleichungen enthalten bei unveränderlicher 
Wandstärke außer den Komponenten des Spannungstensors und seiner Raumableitungen 
nach w je eine der beiden Größen 


1 1 
1 Ar EN B get ß re. 
k > hbiu N RE Re MR DE er 189) 
und 
2 1 
—— Ye ß ne B Fa N TE FE ER; 
u! bi, k u. SER RT (81), 


| welche mithin als Raumtensorkomponeuten mac mg ve 

deutigkeit der Endgleichungen s 

%P und 1® aufzunehmen. Die Auflösung der Gleichungen (80) BE Hu wen 
Ralımen der a ee erg Potenzen der Wan 


Werden noch für die Flächeninvarianten k und I? die Abkürzungen 


BI BERN N ie 4 
B= ne | 


eingeführt, so gehen die Gleichungen (79) über in folgende Form; 


2 a 
x ü BE r Ha he re z g Ha B EEE 33 +b si, +3 87. (85) 
en 
Sl 2 + hlams =. BE, he + em +le(srA| +20: 8rM] 


Zur Veranschaulichung der neu eingeführten Größen sei noch bemerkt, daß t* mit der resul- 
tierenden Tangentialkraft, » mit der resultierenden Normalkraft und m* mit dem Moment der 
äußeren Kräfte in Zusammenhang steht, während n/A einer mittleren Spannung Rechnung trägt, 
welche sich in Richtung senkrecht zur Schalenmittelfläche ausbildet. 


6. Das Gleichgewicht der Schale 


Auf die Einführung der sonst bei Behandlung des Schalenproblems üblichen Spannungs- 
resultanten, wie Längskraft, Biegemoment und Querkraft, wird in vorliegender Abhandlung 
bei Aufstellung der Grundgleichungen bewußt verzichtet, da sich diese Größen nicht unmittel- 
bar auf Komponenten des Spannungstensors oder seiner Raumableitungen zurückführen lassen. ä 
Die Gleichgewichtsbedingungen werden vielmehr — und darin unterscheidet sich die hier ent- } 
wickelte Theorie grundsätzlich von den bisherigen Formulierungen — als Beziehungen zwischen 3 
den Komponenten des Spannungstensors und seiner Raumableitungen nach w gewonnen, wobei 
letztere die Rolle der Biegemomente übernehmen. Auf diese Weise werden eindeutige Flächen- 
tensorgleichungen erhalten, welche davon unabhängig sind, ob der zugrunde gelegte Spannungs- 
tensor kontravariant, gemischt oder kovariant ist. Ferner wird die Gleichgewichtsbetrachtung 
des Schalenelementes im Sinne der Kontinuumsmechanik aus dem Verschwinden der Divergen2 
des Spannungstensors hergeleitet, indem die Integralaussagen 

h 


{ { 
[ Hell wedwe 0. MEET re (86) 
I 


an die Stelle der sonst üblichen, auf dem Begriff des Biege- und Schermomentes fußenden Be- 
dingungen des Kräftegleichgewichtes treten. Dieser rein analytisch deduktive Weg gewähr- 
leistet eine beliebig hohe Approximation an den dreidimensionalen Spannungszustand. 


Die in den Gleichungen (86) auftretende Divergenz des Spannungstensors stellt einen 
Raumvektor dar. Bei der Entwicklung nach der Koordinate w mittels der Taylor-Reihe 


| 
tee] +] Fee 


w=0 =(0 
lassen sich daher die auftretenden partiellen Differentiationen nach w wie bei einem Vektor auf 


= og = se) tw + 

= Sehe I- 

& x Bei der‘ weiteren Pe sind die er, ak A wie Hr ma 
x i 


er und Diuve zu behandeln, so. daß Gl. (67) und (68): zur ır Anwendung Kommen, Unter. Umgehu 
RE ‚einer, Su Zweischenrechnung folgt: ante IERE “ 


Eı Ah= A bs une bP 9 we, Hufe], ie |, u 83] Bro 

DS NT, a bö He ER # se u so sr bs ga man 6 BEN 1 a 5 Au N 2 iR 
y ” RE Ss, ins dB KLEE +90 "5 BEL 5 SE = mr 
are, 2 D san 1. 4902 las... i. ; h RE Hr SD (89) 

ee BR Ra x Ä 
y F 0 F x N « 

ee Lee 

x Yen 2.08 104, —b& Be Rz Hl + 


. den Faktor h°]24 erhält. Wird zunächst das erste Glied mittels der Gleichungen (85) umgeformt, 


so ergibt sich eine Beziehung, we zur Elimination der ES RR KR: Sn JE Ile Baur ER 


Be. mercen kann: 


\ 


msn], vu me N ne . (90). 


ä . In der mit dem Faktor ?/24 behafteten eckigen Klammer treten u. a. die Komponenten srl 


- der Integralbedingung (86) im Falle a = 1 gewonnen wird. Hierbei zeigt sich, daß die mit w 
behaftete eckige Klammer der Gl. (89) für sich gleich Null zu setzen ist. Wenn von den Be- 
ziehungen (85) und (90) Gebrauch gemacht wird und höhere Potenzen der Wandstärke und 

 _—— ‚ihrer Ableitungen, welche über den Rahmen der Potenzentwicklung (89) ande würden, 

vernachlässigt Barden; liefert die Auflösung dieses Ausdruckes nach Du Se 


hl« 


t, 


‚SP 3 | an FR aß 

le 

Nunmehr kann die in Gl. (91) noch nicht berechnete eckige Klammer wesentlich vereinfacht 

werden, da sich die Mehrzahl der auftretenden Größen mittels der Gleichungen (85), (90) und 

(91) eliminieren läßt. Wie die Ausrechnung zeigt, entsprechen die auftretenden Ableitungen 

der Wandstärke h einer Multiplikation des Spannungstensors 8°? mit h, bzw. seiner Raum- 
0 


m’ 
en EN Se Ha + (91). 


3% 


ERINNERN REN EN EITERGENLURIEBSLE TIERE BONRRSELTERN (0, ZUM 


ableitungen RE und S@®l| mit A?/12. Für die weitere Rechnung im Rahmen der durch 
3 3 0 
; 74 =-(Qundn=1 festgelegten Approximation werden daher zweckmäßig die folgenden Ab- 


R- kürzungen eingeführt: 
: ae, gen 
— Ms 


33 0 


h St AR 54° ° gaß 


EN N Rn (92), 


3 


Wird von diesen Abkikkunadn Gebrauch gemacht, so läßt sich Gl. (90) im Rahmen der Ver- 
nachlässigung höherer Potenzen der Wandstärke so umformen, daß die Wandstärke h formal 
nicht mehr in Erscheinung tritt, bis auf einen mit: dem Faktor h?/12 behafteten Restausdruck. 
Wie die nähere Untersuchung zeigt, bei welcher wiederholt von einer in Gl. (61) angegebenen 
Identität Gebrauch gemacht werden kann, finden sich in diesem Ausdruck jedoch nur Größen vor, 
welche in ähnlicher Form schon in den übrigen Bestandteilen der Gleichung enthalten sind. 
Da die einzelnen Terme innerhalb des mit h?/12 behafteten Ausdruckes größenordnungsmäßig 


im Verhältnis :1 kleiner sind als die entsprechenden Terme außerhalb desselben, so kann 


Rare 
120,05 


Bei Anwendung der a (86) auf ea Aucdruck: erhält i im Falle n = 0 das ni RR 
3 Glied (ohne w) den Faktor h, das zweite Glied (mit w) entfällt, während das dritte Glied (mit w2/2) 


() 
auf; diese lassen sich auf Grund einer weiteren Beziehung eliminieren, welche bei Anwendung 


_ “ 
“ 


geführten Tensoren den Schnittkraftresultierenden einer Platte, und zwar liefert T 
"momente angibt. Dagegen besteht keine unmittelbare Identität dieser Größen mit len Schn 


X 


PEST ‚e u 
ri iR ? aA a 14, 5 
He ERW TE Ne 5 


dieser Ausdruck im Rahmen der hier mit » = 0 und n = 1 durchgeführten Approxiı 

vernachlässigbar angesehen werden. Für das ( der tangent a 

‚fläche wirkenden Kräfte ergibt sich dann die a ä- 
u A Zeil, — be Mr), — U. Merl, +63 Men), —bemP — me + = I 


Wie im achten Abschnitt noch erläutert wird, entsprechen die mit den 
der Plattenlängs- und -scherkraft, während M®® den Tensor der 


BE c* : j \Beziehr 
kräften der Schale, welche sich — wie im achten Abschnitt gezeigt wird — aus Bezie 
komplizierterer Art errechnen. ‚Der Vollständigkeit halber sei noch eine weitere Größe ein- 


‚ geführt, welche beim Übergang zur Platte die Plattenquerkraft liefert, bei der Schale jedoch. 
nur einem Bestandteil der Querkraft entspricht: SE | 


er Ba AR E R 5 EURE = % on. 


Mit Bezug auf GI. (85), (91) und (92) ergibt sich, wenn wieder ein mit 12/12 behafteter Rest- 
ausdruck vernachsässigt wird: _ 


x [ a + MePlg To, sea el. DI Reha: e a (95). - Ei Te 
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Bei der Ausrechnung sind die Größen 83] und 2 | wie Komponenten der Tensoren O?#’ 
; 0 3 33 
und. D*#e in Gl.(67) und (68) zu behandeln. Dann folgt: | 
| — 883] — be 8® 4 bp 8 + se seh — b% s=| +52, 8*P 
AO & em 0 0 3 0 & 0 +13 0 


le 
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+ 58 958 — E88 +48] +2jem BE], 4508” 
Bar Br A N a. 13% 22.10 0 |iss f) 
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Die Integralbedingung (86) bewirkt für n = 1 wieder das Verschwinden des in w linearen Gliedes. 
Bei Auflösung nach Ss ergibt sich mit Anwendung der Gl. (85): 
al) 
m* 
ER 2) , 
” * (98). 
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Andererseits tritt bei Anwendung der Integralbedingungen (86) mit n = 0 wieder das erste 
Glied (ohne w) mit dem Faktor } und das dritte mit dem Faktor 8/24 auf, während das zweite 
entfällt. Werden zur weiteren Vereinfachung die Gleichungen (85), (90), (91), (92) und (98) 
herangezogen, so läßt sich wieder eine Beziehung aufstellen, in welcher die Wandstärke % nicht 


‚ mehr vorkommt, wenn von,einem mit dem Faktor A?/12 versehenen Restglied abgesehen wird; 


dieses kann aus denselben Gründen vernachlässigt werden, die bei Aufstellung der tangentialen 
Gleichgewichtsbedingungen (94).maßgeblich waren. Die Bedingung für das Kräftegleichgewicht 
in Richtung senkrecht zur Schale nimmt dann folgende Form an: 


MePlyp ba LP — KM® +n, BEE Pp 0 ne (99). 
Hiermit sind die Aussagen über das Kräftegleichgewicht im Rahmen der füra=0undn =1 
geltenden Approximation erschöpft. Demnach stehen nur drei Gleichgewichtsbedingungen 
zur Bestimmung der sechs Unbekannten L*P und M&# zur Verfügung. Das 'Schalenproblem 
ist daher im allgemeinen statisch unbestimmt und erfordert zu seiner vollständigen Lösung 
die Ableitung der besonderen Bedingungen des Formänderungsvorganges. Hierbei zeigt sich — 


wie im folgenden Abschnitt ersichtlich —, daß sich alle Unbekannten auf die drei Komponenten 
des Verschiebungsvektors der Schalenmittelfläche zurückführen lassen. 


Eingegangen: 11. Juni 1948. 
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Ein Theorem über die 
stationäre Wirbelbewegung kompressibler Flüssigkeiten 
Von Hans Ertel, Berlin, z. Z. Uppsala, und Heilding Köhler, Uppsala 


Es sei & ein Integral des Sysiems der Stromlinien-Differentialgleichungen, ß ein Integral des Systems 
der Wirbellinien-Differentialgleichungen und H die totale Energie pro Masseneinheit im stationären Strö- 
mungsfeld einer autobarotropen kompressiblen Flüssigkeit. Die für Binheits-Differenzen der Maßeinheiten 
gezeichneten Flächensysieme & = const, ß = const und H = const bilden durch ihre gegenseitigen Schnitte 
(&, ß, H)-Einheitszellen, wobei jede (&, ß, H)-Einheitszelle durch den Schnitt eines (&, H)-Einheitssolenoids 
(= Stromröhre) und (ß, H)-Einheitssolenoids (= Wirbelröhre) entsteht. Es wird dann der. Satz bewiesen, 


Wirbeltärke des (ß, H)-Solenoids übereinstimmi mit der Maßzahl der Masse in der durch den Schnitt 
+ dieser Solenoide entstandenen (a, ß, H)-Einheitszelle. 


rential eguations, further H the total energy per unit mass of the steady siream-field of a compressible fluid 
with an autobarotropical relation 0 = 0(P). The surfaces & = const, ß = const, H = const which may 
be drawn at intervalls corresponding to the units of the systems, divide the space into (a, H)-unit solenoids 
and (ß, H)-unit solenoids. T'he (&, H)-unit solenoids, involving the siream-lines, are intersecling the (ß, H)- 

' unit solenoids, involving the voriex-lines, so that (&, ß, H)-unit cells are formed. Then the product of the 
numerical measure of ihe mass-flux through the (&, H)-unit tube and ihe numerical measure of the vortex- 

strength of the (ß, H)-unit tube is found to be equal lo ihe numerical measure of Ihe mass of the (&, ß, H)- 
unit cell. ü 


Soit & un integral du sysiöme d’egqWations differentielles des lignes de courant, ß un ini£gral du systeme 
d’equations differentielles :des lignes de tourbillon et H l’Energie totale par unite de masse dans le ‚champ 
slationnaire de courant d’un liquide autobarptrope compressible. Les surfaces &,='const., ß = const eb 
H = const qui repr£sentent les differences des systömes d’unite forment par leurs coupures mutuelles des cel- 
lules &l&mentaires (&, ß, H). Ohacune de ces cellules prend origine par la coupure mutuelle d’un solenoide 
El&mentaire (&, H) (= tube de courant) et d’un sol&noide Elementaire (ß, H) (= tube de tourbillon). Ensuite 
les auteurs prouvent que le produit forme par la mesure nume£rique du courant de masse & iravers le sole- 
noide (%, H) et la mesure numerique de la puissance du tourbillon du sol&noide (ß, H) est en concordance 
avec. la mesure num£rique de la masse contenue dans la cellule el&mentaire (&, ß, H) formee par la coupure 
des solenoides. ; 


Ilyerp % o6osHayaer uHTerpan CucTeMbI AudpepeHnnansHusx ypasHeHnf NnHnü TORa, 
ß — uHTerpan cnucTemH AuhepeHuuansHEIx ypaBHeunä Anunü Buxpa u H — 06myIio 3Heprmo 
eNeHUIIEI MACCHI B HOJIEe YCTAHOBHBINETOCH TEYCHAA COKHMAEMOH IKUNKOCTH, YAOBHLETBOPAINEH 
YCHOBHW ‚„aBTodaporponnn‘‘ o=.o(p). HosepxHocrn & = const, ß = const m H = const, TOcC- 
TPOeHHbIe MIA HEAIBIX KPATHEIX OT eNeHHM usMmepeHns, o0Pasyl|T CBOHMH NepeceyeHusmu 
eneHnyHBe Kierkn (a, ß, H), mpuyem Kkamıaa u3 HUX OÖpasyercH IepecegeHneM eNeHHYHOTO 
coneHouna (a, H) (= TpyÖku TOKa) U eneHudHoro coneHouna (ß, H) (=Buxpesoä TPpyOkn). 
BATEeM NOKA3BIBACTCA TEOPEeMA, YTO TIPOHSBeNEHHe MEePbI KONHYEeCcTBa SKHAKOCTH, IPOTe- 
KalmımeroO B coneHonne (x H), Ha Mepy Hanpsskenna BuUxXpAa coneHonna (ß, H) pasmo Mmepe 
Macchl B eXeHnyHoH kıerke (x, ß,H), o6pasyemoüä mepeceyeHneM 3TUX CONEHOMNOB. 


I. Dynamische und geometrisch-kinematische Grundlagen 


| Für stationäre Bewegungen einer kompressiblen' Flüssigkeit mit einer für alle Partikel 
| gleichen, Autobarotropie-Relation 


NONE ED eh et De 
(vgl. z.B, V. Bjerknes[l]) hat die Gesamtenergie pro Masseneinheit die Form 


aA dp 
B-Ste+|®......2.2:00...0 


und die auf ein Inertialsystem bezogenen stationären Bewegungsgleichungen (vgl. z.B. 
H. Lamp72] 
RETTET RR N ee ee A) 
bestimmen zusammen mit der Kontinuitätsgleichung 
Rn Cl UNS EN ER EEE a NL) 
und der Autobarotropie-Relation (1) die fünf abhängigen Variablen % = (v,, v,, v,), p und 0, 


während das Potential @ der äußeren Kräfte in (2) eine gegebene Ortsfunktion darstellt. Ver- 
wenden wir allgemeiner ein rotierendes Koordinatensystem (x, y, z) mit dem konstanten Dreh- 


vektor F so ist das System (3) durch 
rm Zuctad He. 2 a2 ee (B) 


daß das Produkt aus der Maßzahl des Massenflusses durch das (a, H)-Solenoid und der Maßzahl der 


Let & be one integral of the stream-line differential equations, and ß one integral.of Ihe vortex-line diffe- 


: en 


a 
4 


Be ersetken! worin P— (6, y; 2) den Ortsyektor bedeutet ind 
das. Potential der Zentrifugalkraft enthält., SE 
Es sei nun «= a (0, 92 ER er u, a Sutemd. RR 


folgt, een sich sofort (vel, 2. B: ‚A, R. Forsyth[3], D. A. Murray) ie Darum 
05 = e„grada xgradH. en 


Sb der Multiplikator e,sich aus der durch Substitution von (10) in (4) ol Bedingung SR Br 


a, EN 5, et 
Gr (m, 92): . | 
als eine Funktion &, —= ex (a, H) bestimmt, so daß nach (10) , 
ES A &, (0, H)grad x grad H I DR I 
euere, wodurch die Kontinuitätsgleichung (4) identisch erfüllt ist. N 
Es sei ferner B= ß (#, y, 2) = const ein Integral des Systems 
| BABES IHNUD FTD En (13) 
der (absoluten) Wirbellinien, welches System der partiellen Differentialgleichung 
u Hm VE RT RE TLA) 
äquivalent ist. Da aus (5) auch 
u won tw, == Re 6. u FREE (15) 
folgt, ergibt sich für den absoluten Wirbelvektor & sofort die Darstellung 
D= en grad B.x.gradH nr a ee (16), 


wobei der Multiplikator eg sich aus der durch die Substitution von (16) in die Identität 


Pag (%) = div rot (v +f x.r)=® 
folgenden Bedingung 


(ce, ß, H) H 
Namur RE N, (17) 


als eine Funktion eg = ep (ß, H) bestimmt; so daß nach (16) die Darstellung 


w='eg.({ß, H) grad ß x grad | 0. u ee (18). 


resultiert. A 


Die für 0% bzw. ib gefundenen Darstellungen (12) bzw. (18) genügen neben den Divergenz- 
bedingungen 


dived)=0, dvW)=0, 
auch den Folgerungen ; 
o% grad =0, %-gradH=0, 


welche aus dem System (5) der hydrodynamischen Bewegungsgleichungen resultieren, so daß 


letzteres, wenn (12) und (18) in dieses System (5) substituiert werden, nur noch eine Beziehung 
liefert, die sich wie folgt ergibt: 


a nee re Me 


——n me 
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Zunächst ist 
| &, ep (grad & X grad H)x (grad xgradM)=ogradH. ...... (19) 
oder: \ 
& ep { (grad & x grad H) grad H) grad ß— ((grad x grad H) grad ß) grad HY = ograd H, 
"und wegen ’ 


(grad x grad H)grad A =0, 
sowie | 


2% (0, ß, H) 


‚(grad a x grad H) grad = — FRE 


sind dann also die hydrodynamischen Bewegungsgleichungen (5) mit (12) und (18) erfüllt, 
wenn die Beziehung - 


| ; (% Y; 2) 2 


besteht, welche wir, die absoluten Beträge einführend (oist immer positiv), 


Io (6 Bil os, De EL, ee 


schreiben. Diese Gleichung (20) ist der analytische Ausdruck für die funktionale Unabhängig- 
keit der drei Funktionen «, ß und H, welche Tatsache in geometrischer Interpretation bedeutet, 
daß die .äquiskalaren Flächensysteme € 


8.12% 2) const 
E: Piz 2) const‘; 
H (x,y,2) = const, 
sich gegenseitig schneiden müssen. 


Denken wir uns diese äquiskalaren Flächensysteme gezeichnet für 


GA NE EV ESCONSE en Te en EN (21), 
N RE RE a TEE BEN MER Der (22), 
HEHE Foster Ba (23), 


worin Ey Es, Ep die Maßeinheiten der Größen a, ß, H bedeuten, also z.B.: Ey = 
1 erg/gr, und n,, N8, Na =1,2,3,4...,' so wird das Strömungsfeld durch die «-Flächen in 
(«)-Einheitslamellen, durch die ß-Flächen in (ß)-Einheitslamellen und durch die H-Flächen in 
(H)-Einheitslamellen zerlegt. Durch die Schnitte der «- und H-Flächen entstehen die (a, H)- 
Einheits-Solenoide, die zugleich Stromröhren gemäß Gleichung (12) darstellen, während durch 
die Schnitte der ß- und H-Flächen die (ß, Z)-Einheits-Solenoide entstehen, die gemäß Gleichung 
(18) Wirbelröhren sind. Die Querschnitte der (x, H)-Solenoide senkrecht zu den Solenoid- 
achsen sind durch 


OH Dy— Dr Ea}]| grade X. grad HN DR RR 
die der (ß, H)-Solenoide- durch 
AN RE VE RE EL St 2 ER Y5)) 


gegeben. Schließlich wird durch alle drei Flächensysteme (21, 22, 23) das Strömungsfeld in 
(a, ß, H)-Einheits-Zellen aufgeteilt, wobei jede (a, ß, H)-Zelle auch durch den Schnitt eines 
(«, H)-Solenoids mit einem (ß, H)-Solenoid entstanden gedacht werden kann; Das, Volumen 
einer (&, ß, H)-Einheits-Zelle ist durch 


9 (0, ß, H) 


0 (x, Y,2) > 


V(&B, H) = Ba: Ep: En] 


bestimmt (vgl. hierzu z.B.: H. Ertel [5], [6]; für Potential-Strömungsfelder hat bereits 
M. Lagally [7] ähnliche Betrachtungen durchgeführt). 


E 
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VE heorem her MR Bekietung 


an die Masse dr (im ehe | 
sich aus (12) und (24): 


röhre mit konstantem Wirbelmoment darstellt. 


en en ET rt 
See ee pl : er aan und de M 
KIEERIEEE bildete: 0,8, H )-Einhe we Bu 


Aus (20) und (26) folgt zunächst durch Elimination der Funktionalde termin 
ea RI Berka in > 


00 DI= 8181016, M) (Er En) = FM) (Be En) en. (AR 

ER F(&, H)= 0|%|Q(& H) den Massenfluß durch das DSB verstanden, und * 

analog folgt aus (18) und (25): Wr 

[es (ß» = BO [En = W (BE) IE) Sr 

worin W (ß, H) = |%| @ (ß,-H) das Wirbelmoment des (ß, H)- Br ua darsiellt Die Subst A: 

tution von (28) und (29) in (27) ergibt dann: + E 
Fa, A): Wi, = MB, H)- Ep... .\. 


| * Rs 
Wenn nun weiter F (a, H), W(ß, H),und M (c, ß, H) die Maßzahlen der ee den Ei 
Einheiten E7 (z. Br —= 1gr/sec), Ey (z.B. = 1cm?/sec) und Ey (z.B. 1gr) ausdrückbaren 
Größen F (a, H), W (ß, H) und M (a, ß, H) bedeuten, so ist also gemäß der Definition einer 
Su Größe als Produkt von Maßzahl und Maßeinheit: 


besteht, und die Substitution von (31), (32), (33) in (30) führt dann mit Rücksicht auf (34) 
zu folgender Beziehung zwischen den Maßzahlen von Massenfluß durch ein («, H)-Einheits- 
Solenoid, Wirbelmoment eines (f, Z)-Einheits-Solenoids und Masse der durch den Schnitt 
dieser Solenoide gebildeten («, ß, H) es Zelle: 


ram tom-ham) 1 ER +: | 


womit folgendes Theorem bewiesen ist: 


F («, H) I H): Er BIER DI RFR I» EN (31), 
WIR. el BD N 32,208 
M (68, = Mila BED Bir a he (33), 
während zwischen den Maßeinheiten Ey, Ey, Ey und Ey (= lergj/gr) die Relation, 
Bu Dr er SB ET (2 | 


Theorem über stationäre Wirbelströmungen: 


Das Produkt aus der Maßzahl des Massenflusses durch ein 
(die Stromlinien enthaltendes) (, H)-Einheits-Solenoid und der Maßzahl 
des Wirbelmomentes eines (die Wirbellinien enthaltenden) (ß, H)-Einheits- 
Solenoids ist gleich der Maßzahl der Masse, welche in der 


durch den Schnitt dieser Solenoide entstandenen (&, 8, H)-Ein- 
heits-Zelle enthalten ist: 


Gorollar 1: Für verschiedene (a, ß, Z)-Einheits-Zellen desselben (x, H)-Einheits, 
Solenoids sind die Maßzahlen der Wirbelmomente den Maßzahlen der Masse der (a, ß, H)- 
Einheits-Zellen proportional. 


Beweis: F (&, H) = const längs desselben (a, H)-Einheits-Solenoids, das eine Strom- 
röhre mit konstantem Massenfluß darstellt. 


Gorollar 2: Für verschiedene (a, ß, H)-Einheits-Zellen desselben (ß, H)-Einheits- 
Solenoids sind die Maßzahlen des Massenflusses den Maßzahlen der Masse der (a, ß, H)- -Ein- 
heits-Zellen proportional. 


Beweis: W (, H) = const längs desselben (ß, H)- )-Einheitssolenoids, das eine Wirbel- 


‚ 
N er 
nenn nn nn 2 


nn + 


! 


40808 eigen =1- ‚Em aber ‚wir unterdrücken | im Pigenden 4 der . Kürze. halber die Mat 
en so resultiert: ? EIRSE KERN | 4 2 
A er RR 
ER ER a n l 6,% = — 6 ı EN Er 
ön: il, 12 | - sin (8 D) ER om er 


Baar n zerer 
I f ur 1 A « x i 


Meat Heim a H. 3: am HE 2.2.0.8: 244). Nü un kann das Volumen einer baratierkpipedf: .; 


& schen (&, ß, H) Binkäits- Felle, mit Hilfe ‘der Kantenlängen sa Ka EZ ee \ 
R Sen m Br H.Schnitikante) in eat Weise dargestellt ER RE 2 


ER ag ER RN} sn nd m): su = Via P, A) ee Rd 
Banner. en = VB, m R ee... 
5 EEE Sam? sBn* sin (®, ») dIn=V 1% B, A) m TER RER, EN 


Multiplikation yon (36) mit. (37) ind (38) ergibt nun: 
#1 Q1 a, H MOB. m: Sr Sp Sin (, Don pr (aß, N 


) und nach Multiplikation init der Dichte oe und Division durch ( (39) folgt: 


Pe: elil0“ MOM MuEM-ermhM 
Er oder ; SR 


f 


FM) W(ED)=M (aß, m Be a) 


welche Gleichung | mit. ( (30) bzw. (35) bis auf die-hier nicht abgespaltenen Maßeinheiten iden- 
SCH Ist... o.- 
Der Vorteil des Theorems (35) Besenäber der Gleichung (36) liegt darin, daß dasselbe nur 
nn miteinander verknüpft, während die Gleichung (36) vier Größen (ön, |%], j%], 
sin (2, %)) enthält, und man erkennt, daß diese Reduktion der vier Größen auf drei durch die 
NE des geometrischen Begriffs der («, ß, H)-Einheits-Zelle ermöglicht wurde. Es ist <A 
ohne weiteres ersichtlich, daß obiges Theorem auch für die stationären DIEBE DEIUNE ER er. 
an unpesalbles Flüssekeifen gilt. ER u Me 
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‘ dem man abschnittsweise ‚‚durch‘ je N = 2n +1 aufeinanderfolgende Punkte &;_„, -.- Zp4n 


5 ki a ; 7 | i , Ev 3; 4 
 Mehrfaches Ausgleichen einer fehlerhaft 
\ 0... Von K, Stange in Bad Schsa 


"Das mehrfache Ausgleichen einer Punktreihe xx (tx) nit EEE. 
Kal en Ole wird bereits nach ) nahezu ie U e mit jedem). 
' Schritt mehr und mehr wachsende Korrelation zwischen den Fehlern u” der geglätteten Werte x. Außerdem 
a ER ee : a ae tr er, y ee igerwren.. 
ie einfachsten Fälle linearen 1 jet. A 
De iin der Hintup einer von vornherein vorhandenen Korr der Ausgangsfehler u ar 
fortpflanzung, untersucht, ER; x . 
The iterative graduation of a set of points wx (tk) afflicted with errors, by add of power series 0j a limited 

degree proves ineffective almost after the very first steps. The cause is Ihe correlation between the errors u" 
of the smoothed increasing more and more with each step. Moreover, the primary function = Ba 
“jalsified” at each step by am ‚amount. The conditions are computed in the most 


FW 


w 


cases of linear and quadratic . Besides, it is investigated, how a correlation of the origimal errors u 
existing from the first time is influeneing the propagation of errors. 3 nt KIT 


de d’une compensation reiteree d’une serie de x (tx) chargee d’erreurs ü Vaide dune 
RE ie and de degre limite devient 1 Tape "a Cu na a ro ee in 
reurs u‘” des valeurs des qui s’aceroit de . En outre, la fonction initiale x (t) est de m&me 
faule! L’auteur donne le calcu ee d’une compensation lintaire et carre. De 
plus, il montre V’influence d’une correlation des erreurs initiales u sur la propagation des erreurs. 

Merox nocHeAoBATeNbHBIX TONPABOR HEheKTHOTO pana TOyerK 2 (l£) IPH IIOMOIIH CTEeTIEH- 
HEIX PAAOB OTPAHUYeHHOH CTeIIeHH CTAHOBUTCH HEMNeÄCTBHTENILHEIM Ye HOCIHIe HECKONBEUX 
mmaTOB. ITO OÖBACHAETCA PacTyımek C KAKABIM IIATOM Koppesmmmeili mekay OIIHÖRaMH (9) 
CTIIAKeHHEIX 3HAyeHHÄ . Kpome Toro, B Toä-xe mepe HUCKA:KAETCH © KASEIEIM ITATOM H 
HcXxonHas hdyHRuuAa © (6). ITu YCHOBHA YUIEeHHO HCCHEAYIWOTCH B IIPOCTEeHIIAX CHYYaAX AHHeii- 
HOTO H KBANPATHYHOTO MeToya MoImpaBor. Kpome Toro uecaenyeTca BIHAHHe umewmelich € 
CaMOTO HAYAIa KOPPEeNIANNUH Meky HCXONHEIMM OIINHÖKAMH u Ha PAcpoCcTpaHeHue OmMOOK 
B Ipomecce HONPABOR. . 


.1. Einführung. ‚Glättet man eine Reihe von Meßpunkten ... 2g—1, 7 Zi, - ++, IN- 


die „beste“ ganze rationale Funktion vom Grade A legt, so ist der ausgeglichene Wert z, an 
der Stelle }, 


4. 
»-= > > Zr I Pr Re (4): 
> vn v 

}, sind die Gewichte, mit denen die fehlerbehafteten Ausgangswerte x,., in die Berechnung des 
‚„‚Mittelwertes“ eingehen. (Sofern nicht ausdrücklich etwas anderes bemerkt wird, soll der 
Summationsbuchstabe » stets alle ganzen Zahlen von v»= —n bis »= +n durchlaufen.) Ei 
Falls die so gewonnenen Punkte x; noch zu stark streuen, liegt es nahe, das Verfahren zu wieder- 
holen und aus x; die zweimal geglätteten Werte 4 


zu bilden usw. Bei großen Punktzahlen N und höherer Ordnung A der Ausgleichfunktion wird R 
die mehrfache Durchführung allerdings wegen des großen Rechenaufwandes „praktisch un- 
möglich“. Es gibt jedoch zwei Ausnahmen, nämlich die Ausgleichfunktionen ersten Grades 
(?£=1) ohne Rücksicht auf die Punktzahl N und die Ausgleichfunktionen zweiten Grades 

(R= 2) für N=5, wo der geglättete Wert %, sich auf die einfache Form 


3 
| 35 Ak RZ SER ea en A 
bringen läßt. In diesen beiden Fällen scheitert das mehrfache Glätten jedenfalls nicht am 
Rechenaufwand und wird zur Herabsetzung der Streuung der x-Werte auch tatsächlich emp- 
fohlen ),. Da die Gesetze der Fehlerfortpflanzung beim mehrfachen Glätten in Abhängigkeit 
von der Schrittzahl a offenbar noch nicht untersucht worden sind, solles im folgenden geschehen. 
Dabei stehen grundsätzlich zwei Wege offen. Entweder drückt man auch die o-mal ausge- 
glichenen Werte «% durch die Ausgangswerte ..., %, ... aus. Dann darf man das Fehler- 
fortpflanzungsgesetz in der üblichen Form benutzen, da die Fehler der Meßwerte x, als un- 
abhängig voneinander angesehen werden. Oder man berücksichtigt die zwischen benachbarten 
Fehlern der geglätteten Werte bestehende Korrelation. Dann hat man das Fehlerfortpflanzungs- 


x, = % — 


!) Fr. A. Willers: Methoden der praktischen Analysis. Berlin und Leipzig 1928. S. 252 und 253. 
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ey ist leicht zu übersehen, daß der Summationsbuchstabe o nach o Schritten vne=—n 
bise= + no läuft. Es ist jedoch für das folgende zweckmäßiger, oe unabhängig von der Sch ! 


Mu a Ber Cr } Ken BEE aA Sort 2 $ ’ 


Bi I Die che ., a deneh, die Bat Fe, in "die ne Mitte 
| e” eingehen, werden durch zwei Zeiger o und o gekennzeichnet, von denen sich g auf die 
des x-Werts und o auf die, Schrittzahl des Glättens bezieht. ‚Faßt man beim Ken 


r 


zahl d stets von — & bis. + © laufen zu ‚lassen und diese Festsetzung durch N AR NEE 


Lee 


zu ergänzen. Für 6 = 1 ist en RT ie — m = e< = N ung au nz B Für = is 


kann sind. Es ist nach (3) RN ya Hr | 


Anschaulich geht das Bildungsgesetz für die Gewichte aus Bild 1 hervor. Das Gewicht Pa, für 


h Er paT Et je 
Mi; re = a are RAR 
Rt Y e=—on z Bar ö "Fa 4 


y an a 


N=2n +1 Punkte ‚zusammen, d.h. wählt man je n symmetrisch zu &, ‚liegende x-Werte a 


% 


N Pe für en WE m 
0) 


a 


dürfen wir also e durch » ersetzen. 


Ir berechnen nun die Gewichte en, 


‚ wenn u des yorausgehenden Schrittes u ben 
Ra) = Z% RCB < . ? ee EL 1e s ; a 


und Bach (A ‚) mit der oben getroffenen ER über des Ehteich. von 6, “ | 
. ‚(o) h 


; Tk-Ly Zah Bine = he X:-+o ER a ER a Bd. 


Setzt man den letzten Ausdruck in 5 ı) ein, so Bone VE Ren 
EI=LH, Sera - | BL: 

oder nach Umordnung der Summe, | \ | 
ar 2(&% he) herren a . (6,). 


0 


# 


Andererseits ist macht (4). er | 
£ zer _ set Veto ee a ae EEE eh We (6): 
j . e 
Durch Vergleich von (6,) und (6,) gewinnt man die gesuchte Rekursionsformel 


4 BR HLN Mu) = Shhens at, 


irgendeine Stelle x;+, ist nach (o +1) Schritten gleich dem IRRE) zu (1) gebildeten] 
„„Mittelwert‘‘ der (2n +1) symmetrisch zu 
2 liegenden Gewichte A „,... 20 „des 
vorausgehenden Schrittes o.. Man glättet 
also gewissermaßen von Schritt zu Schritt 
die Gewichte. Da die Ausgangswerte Ay 
durch eine gerade Funktion von o gegeben 
sind, AD =A,—=4_,= A), so werden nach 
Bildi auch sämtliche übrigen Gewichte 
gerade Funktionen von o Bild1. Zur Entstehung der Gewichte at) 
An — A) RR beim mehrfachen Glätten. 


(6) 
A 


(0) (0) 
ER. 


9 Kl 


Symmetrielinie für die Gewichte A'9? 


ir 
T 


gt 


SER en 
punkt x; ba Teilbereichs bezogenen Nomente 


der Gewichte 4 maßgebend. Mi” entspricht dem statischen Moment und ar il r 
moment der Gewichte. | Pia 


falls « ungerade ist. Weiter wird 
' Setzt mn o—r= es = @ +2ve 7 »® und Be die letzte Doppelsumme um, so kommt 


| oder mit (9;) 


h: 3. Nerehnung dr Momente PR re 
r die „Fälschung“ f, der 


W=LEM, nam >0. “u ur A > 


a Zunächst ist für «= 0 gemäß (7) ur Pe: et 
sr = anmz -2(2 RS, )= ze EnM- un gn 


e 
Da nun My RN ist, so wird Pi 
y v k ° u” % 
Mm -2 20 — = { ” ” h ” E E ” * ” L .22W 3 ” - Kir e - 


Die Summe aller Gewichte A) hat nahen von der Schrittzahl o den Wert 1. ' 
Da A) gemäß (8) eine gerade Funktion von og ist, so verschwinden sämtliche Momente 
ungerader "Ordnung, also e 


S 


Ms Mezeu = ZeLhh,. 


MEN S [2 2. PETE »zar) 2 
Re Yv e € 1:3 


ms) 3%, (nor 2»M’+®Mt. 


Daraus wird schließlich wegen (9,) und (9;) 


N 
M—=MP+M) ode, M=oeM)........ 9a) 

Eine ähnliche Rechnung liefert für M"” die Rekursionsformel 1 
| MIZMT LEMON MO T MEN EN Eu ()- 


Im Hinblick auf die eingangs angedeutete Beschränkung der Ordnung der Ausgleich- 
funktion auf die Werte R=1 und R= 2 .können wir die Berechnung der Momente mit M, 
abbrechen. 


4 

4. Mehrfaches Ausgleichen mit Hilfe einer Geradenschar. Bei schwacher Krümmung j 

der Funktion x (t) glättet man sie abschnittsweise mit Hilfe der ‚‚besten‘‘ Geraden „durch“ | 
je N=2n +1 Punkte. Die Gewichte A, sind dann unabhängig von » N | 
a: | 

Er es er Ve (10), | 

so daß bei mehrfacher Anwendung des Verfahrens gemäß (7) | 
Re > a | 

I Rye A ns Mr UN RE (11) 1 

wird. Man findet also das Gewicht AY*” für den Wert %;+. In der Zeile (oa +1), abgesehen 1 
von dem Faktor 1/N, als Summe der N=2n +1 symmetrisch zu %; +. liegenden Gewichte . 4 


für N=3 gebildeten ‚Dreieck der Einflußzahlen“, von dem selbstverständlich wegen der 
Symmetrie die rechte ‚Hälfte‘ genügt. 


’) K. Stange: Die zweckmäßige Auswertung von punktweise aufgenommenen Zeit-Weglinien. 


2, der vorausgehenden Zeile o. Das kommt anschaulich zum Ausdruck in dem folgenden | 
Ing. Archiv XVI (1948), 8.391 u. 392. | 
j 


; 2 = 2 geht Are; Gewichtsverteilung für us 02 
ersten Schritte aus Bild2 hervor. Die „Gewichts- 5 
BR kurve“ ‘verbreitert sich mit wachsendem o mehr 

vE und mehr, wobei jedoch det, Flächeninhalt zwischen 

=3,Klırve‘”, und g-Achse gemäß‘ (9,) Bande gleich. B: 

20; we - SE 

Per “alle; alle Were mit der eh Genauigkeit 
gegeben sind, so wird der nach o Schritten noch vor-ı, 


Bi . handene Fehler“ Rz Band : 
go £ RR ad. 


Aue wobei go der miitlere ne der Rev aakfe x 
Bist... Das Fehlerverhältnis ist in Bild 3 über 
E90 dargestellt. Wir haben uns dabei auf die einfachsten 
= Fälle N=3 und N = 5 beschränkt. N = 3 entspricht: 
einem oft benutzten zeichnerischen Verfahren ®) und 
00. N=5 ermöglicht einen Vergleich mit dem Fehler- n > % 


= — abfall beim quadratischen Glätten gemäß (2), woeben- zia2. Die auf die Ausgangswerte ae 
‚falls fünf Punkte zur Bestimmung der Ausgleich- zogenen Gewichte A(® beim mehrfachen Glätten 


E ”; - parabel erforderlich sind. x mit Hilfe der besten Geraden durch je 5 aufein- 
B-: ; anderfolgende Punkte in Abhängigkeit vonder 
.Die infolge des begrenzten mathematischen Schrittzahl o. 5 


 ——  Aufwands beim Glätten in Kauf zu nehmende Fäl- 
F schung fm der Funktionswerte wird nach co Schritten in erster Näherung *) 


Z 2) — Kg RN TEILEN. IE (13,). 


2! 


Bild 3. Abnahme des mittleren Fehlers beim ‚Glätten mit Hilfe der besten Geraden (R = 1) bzw. Parabel 
(R = 2) durch je N aufeinanderfolgende Punkte bei mehrfacher Anwendung des Verfahrens. 


%;, und a sind die Sähren von Meßfehlern freien Werte der Funktion und ihrer zweiten Ab- 
‚ leitung an der Stelle i,, r ist der unveränderliche ‚‚zeitliche“ Abstand der Meßpunkte und x) 


nl Willers: a. a. 0. 8.252, Bild 107. 
4) Stange: a.a.0. 8.392 Gl, (23,) mit o an Stelle von v. 


beeinflußt wird. Insbesondere wird für N = 3 bzw. n— 1 


Ai RN Pe TREE A n 
Be Bei jedem Schritt fälscht man ie Funktionswerte demnach um tens 


(tl) eine lineare Funktion, für welche die zweite 
man das Ausgleichverfahren beliebig oft wiederholen eig | 


3 R ER er. ER © EN 5% [u 5 
Smd für N 5. baw. wide 2 | Ba? 
pe (2) 2 OF ze = az (1) BL Er ii De te . . 2a), > Ai: Er 


d.h. bei Verwendung von je fünf Punkten fälscht man die Ehen; x (f) mit jedem Schritt 
dreimal so stark wie beim Ausgleich mit nur je drei Punkte 


5. Mehrfaches Ausgleichen mit Hilfe einer Parabelschar. Bei stärkerer 
der Kurve & (£) wählt man zum Ausgleichen Ersatzfünktionen zweiten Grades. Sie liefern beim 
' ersten Schritt für N = 5 den geglätteten Wert (2), den man mit Hilfe der Gewichte ), in e= 
für die Bund geeigneteren Form 


H=& 4, mit =, (175%) ne: (14) 0 
S, j 


darstellen kann °). Aus (7) berechnet man die Gewichte’A/” und die für die Fortpflanzung der 


Meßfehler fl maßgebende Summe ihrer Quadrate. Man findet für die ersten Schritte den in 
Bild 3 dargestellten Fehlerverlauf. 


Zur Berechnung der Fälschung f}? beachten wir, daß bei einer Ausgleichfunktion zweiten 
Grades das Moment M5 $ verschwindet 9. Die Rekursionsformel (9,) vereinfacht sich dann zu 
MY=- Mm" LM” wundliefrt MO=oM" ...... (15,).. 
Aus (9,) folgt 
MEN NER (15,. 


Bei einer Ausgleifhfunktion zweiten Grades verschwinden die Momente zweiter Ordnung un- 
abhängig von der Schrittzahl co. Der mathematische Fehler wird demnach in erster Näherung 


” 1 
Mao 2.8 5 Wir MO EN LTE De (16,). i 
Mit (15,) und MV -_ Be — re — sun +1)(n+n—2) wird { 
y 30 lr. 2 f 
ae A USE A "er (n r)* (q + fi +45) en N (16,). { 


Mit jedem Schritt fälscht man demnach die Funktion um den gleichen Betrag f{. Funktionen 
bis zur dritten Ordnung einschließlich, deren vierte Ableitung verschwindet, bleiben ungefälscht, 


auch wenn man das Ausgleichverlahren beliebig oft wiederholt. Insbesondere wird für N=5 
bzwin ='2 


d 3 
ee ELLE 


und für‘ N = 7 bzw. n=3 
a) BE 2)... RS ET Sr (16,). 


°) Willers: a.a.0. 8.253. Gl.(11) und (14). 
°) Btange: 3.0, 8.387, Gl.(11,) und 8. 394, Gl. (28). 


UT 


r Br selmaiker De 
n da Verfahren jedoch wie r 
aı rbundene ‚Fälschung der 'Funktio a so ge sche FR 
tte rnoch. eine geringe und schließlich fast gar keine Verbesserung der Stre 
Infolge der „langsamen. Bea des. ‚Verfahrens ist, es ‚nicht. an 


ortzusetzen, bis.die Verbesserungen - = > A unmerklich werden. Über! eo. 


zugehen, ist im allgemeinen zwecklan® Die aus Bid: 3 ersichtlichen Eigebnissä‘ werden düre 1 
“ die Erfahrungen des a ig Be ee an ar in vol 
De sn bestätigt. 
Mit‘ Hilfe von (13,) Sad ‚(16 5 ae man imstairde, den NR o Schritten ER 
ar Sr ea Fehler fi völlig zu beseitigen, wenn die wahren meßfehlerfreien Werte &,(f) 
bzw.'&, (£) der zweiten bzw. vierten Ableitung der Funktion bekannt wären. Das ist natürlich 
So bei Erfahrungsfunktionen x (£) nicht der Fall (sonst wäre das Ausgleichen auch unnötig). Der 
Wert der Gleichungen besteht also zunächst darin, daß man die Größenordnung der Fälschung 
richtig abschätzen kann. Darüber hinaus gelingt es aber, die Fälschung Fin‘ wenigstens. Be 
M großen ganzen‘ herauszubringen, wenn sich die Ableitungen & (t) und '%& (t) mit einiger Sicher- 
heit bilden lassen, was tatsächlich oft der Fall ist. Man ermittelt dann nach Abschluß des Aus- 
gleichverfahrens aus der geglätteten und gefälschten Funktion x (t) die erforderlichen Ab- 
‚leitungen, glättet sie notfalls noch einmal und setzt sie als „Ersatzwer “ für die unbekannten 
ar genauen Werte &, (f) und &, (f) in die Formeln (13,) bzw. (16,) ein. Verbessert man. die aus- 
Br: Behlicheren Werte x‘ (t) um die so gefundenen Beträge f(®, so gewinnt man eine Funktion 7 
Be} ”.(), deren Streuung mit der Streuung der geglätteten Funktion x (t) übereinstimmt (also - 
eher als die der Ausgangsfunktion « (f) ist), die aber andererseits im Gegensatz zu 9 (t) 
nur noch -um ganz unwesentliche Beträge vom wahren Verlauf &, (t) abweicht. Grundsätzlich ' a8 
kann man die Fälschung mit Hilfe eines Iterationsverfahrens noch etwas besser beseitigen, indem 


z man von «2 (t) erneut die Ableitungen bildet, diese in (13,) bzw. (16,) einsetzt und damit = (t) 2 
RR verbessert usw. Im allgemeinen wird sich dieser Aufwand aber kaum lohnen. Kr 
ee Beim linearen Ausgleichen läßtsich für N = 3 der geglättete Wert z{+" — ze, +. 
+ +21] auch in der Form ’ 
a RL 
2 Ra — 20) re RE ER NE: (17,) 
E ae 
A darstellen, welche der Beziehung (2) beim quadratischen Glätten RER 
E on 5 B 3 t i ’ r 2 = 
E’ { wet & — 35 AN: US RR st (17,) | ; 
entspricht. 4x und A4” sind die aus den Werten x, gebildeten Differenzen zweiter w e 
bzw. vierter Ordnung ‘an der Stelle i,. Andererseits gilt für die entsprechenden fehlerfreien Ba 
% 


I 


Differenzen 
KR 2 Fa 4 ae 
HE Sand FA 


'so daß man (13,) und (16,) auch in der Gestalt ‚ 
| DIT EN bzw NOS SER ee SEES 
m 3 [3 . m \% 35 * ? 1; 2 


schreiben kann. Vergleicht man nun (18,) mit (17,) und (18,) mit (17,), so sieht man, daß die 
zur Herabsetzung der Streuung anzubringenden Verbesserungsglieder die gleiche Größenordnung 
wie die in Kauf zu nehmenden mathematischen Fehler fi haben. Das ist auch anschaulich 
einleuchtend. Denkt man sich nämlich das Verfahren so oft durchgeführt, d.h. o so groß, daß 
die Funktionswerte «(9 (£) „‚praktisch fehlerfrei‘‘ geworden sind, so würde die Fortsetzung des 
„Glättens“ nur noch eine Fälschung der Funktion hervorrufen, die aber mit Hilfe der Glei- 
chungen (18, 9) wieder vollständig rückgängig gemacht werden könnte. 


’ 


al 0 SENDE Lina Gala rag, 32) ABER N le lan Da md nn 2 


ee aus der IR Aus er lee ahatch in diesem Ile 1 


Rechnen wir nun nicht „‚sprunghaft‘“,sondern „stetig“, so genügt die Genichtstunktion (6; 


Ne des zweiten Schrittes, a 


KRK=+b+22...) voneinander unabhängig sind. 


el it ee der ; 


>R „'n 


AN x a Tr +M HaRa=at Eu 2 
oder mit Hort ER | EEE ; nA 
/ an | ziern _ ya 1 224% RE RL 
Be ER aha Re Pace EI ER Dar TE 2 N 
do- 3: de uud 


der partiellen Differentialgleichung 


aa 1984 Be | 
a0? Be | 


die in der Theorie der Wärmeleitung 4 eine Rolle spielt und deren allgemeine Lösung Be 
ist ). Wählt man als Ausgangswerte A (0;0) die durch eine ange dargestellte 


1-58) für jel<3 und 180) 20. fie, Mas SE. 


so stimmen die stetig und sprungweise berechneten Gewichte praktisch überein. Man it so 52 
imstande, die Gewichtsverteilung für irgendein o = co’ +2 zu berechnen, ohne die Zwischn- 
stufen alle durchlaufen zu müssen. Da jedoch in unserem Falle im allgemeinen nur die ersten 
Schritte von Bedeutung sind, so soll der angedeutete ig nicht weiter verfolgt 
werden. 


79. Die Korrelsiionsfunktion Te Wie im Abschnitt 1 bereits betont wurde, ist das 
Gesetz für die Fortpflanzung des physikalischen Fehlers beim mehrfachen Ausgleichen nicht 
leicht zu übersehen, da benachbarte Fehler schon nach dem ersten Schritt nicht mehr unab- 
hängig voneinander sind. Wir bezeichnen im folgenden mit zw, die Fehler der Ausgangswerte #,; 
und mit x, die der einmal geglätteten Werte x,. Berechnet werden soll das Korrelationsmaß 
(im folgenden kurz K-Maß genannt) 7,, welches zwischen ‚‚benachbarten“ Fehlern 7, und fr+, 
im Abstand o(e=0;+1;+2;...) besteht, wenn ‘die Ausgangsfehler #, und 444, 


Allgemein wird das zwischen zwei Reihen £&; und 2 bestehende K-Maß r - durch die Glei- 


chung 
& EN 
3 RR (20) 
ZI n: 4 
erklärt ®). In unserem Falle wird wegen a ms— Si ü, der gesuchte Wert 7, (bei festem } 
e=0;+1;42;...) A 
& Üp are e; ; 
Be tn RE a ER (21) { 
; 
Nun ist auf Grund von (1,)' » 
1 = - Aylipyy Und fire = - Iylintetv nn nn ee (22). 


Setzt man diese Ausdrücke in (21) ein und ordnet die are dreifache Summe um, so 


erhält man Y 
fi 
Ber = Ayhys '(@ Mx-+v TER) . 


he ILhhr (Hot) R R & 

Wegen der Unabhängigkeit der Ausgangsstreuungen ist im Zähler bei festen Werten von » und »’ | 
0 für o+v’ #9 

2, Mick Pirachv” == B u für o-+ en a Arte SE (23,) y 


.) B. Riemann: Partielle Differentialgleichungen. Braunschweig 1938, 8. 110 und folgende, 


1940, AN a Hugershoff: Ausgleichsrechnung, Kollektivmaßlehre und Korrelätionsrechnung. Berlin 
) 


, u 


Be NR . Hr 3a h von. 


Oel SAN 5 
ei 


4 für at u wird schließlich 


| a 
Stets ORDER Bi - 


Al md 7-0 für Be rn OR 


"Wir wollen 7, ‚‚die durch das Ausgleichen erzeugte Korrelation“ nennen. 


Bildet man die Summe aller K-Maße und beachtet, daß man o wegen (a: ) bei der Sum-’ 


_ mierung von — oo bis “ co laufen lassen darf, so wird 
E +2n 


N 1, ER a Er 
en Den, en hrs 2 


"Nun ist- bei festem y 


© + oo 
2 Ayte ee 
: . 0e=— » K=—o =—n 
Also erhält man 
ö 5 +2n 1 
> un = 7 RL ie ee ee . (25). 
, e=—2n N) 


Die ne aller K-Maße ist gleich dem Kehrwert des Gewichts A, in der Mitte des Bereiches 
der zusammengefaßten Meßpunkte. 


Beim linearen Glätten mit Hilfe von je N=2n +1 Punkten erhält man aus (24,) wegen 


1, = konst = 2 | 
SE Se ee) 
v=—n-+o N N 


Die K-Funktion hat also „‚dachförmigen‘“ Verlauf. 
Beim radrabischen- Glätten mit Hilfe von je 5 Punkten folgt aus (24,) ausführlich auf- 
geschrieben ®) 


ee EEE ER u nt 
Ar RE NEN) = h, EN) 0,565 
0 
RE RS RAN. ro 
y I (27). 
he eh) ER N 2 oa 
A) 2 
| 1.43 
u Mai) Sk = -- 0,015 
e | ur a0 


Es besteht also nur zwischen unmittelbar benachbarten Werten ji, und ft. Korrelation, 
während die Fehler mit größerem Abstand trotz des Glättens nahezu unabhängig voneinander 
bleiben. 


) K. Stange: Das Bildungsgesetz für die Fehlerformeln beim Ausgleichen von WERreihen mit Hilfe 
ganzer rationaler Funktionen wachsender Ordnung (nicht veröffentlicht). 
10) Die Zahlenwerte 0,565 usw. ur man mit (14). 


Damit ist Has K-Maß 7, ieln benachbarten Berl ii, a ME in NDR Weise dureh 2 
die bei der Mittelbildung an Gewichte A, EIERN! Mit 2) an }, und 4 =0 RER 


= 1, he; ei<2r a en 


im Vornusgehöäden A > Sacl 
‚auf die Aufgabe, das Gesetz für die Fehlei anzut 
berechnen, wenn die Ausgangsfehler wy nicht mehr uı 
ige sind. Im folgenden sei also das K-Maß fe‘ 
ir und Hirte als Funktion von o bekannt h he > | 
in Mit ERRSEHN | ve z = Fi REN 
== Vesz n 7 —=i1 i one en al Sc 
n= Pr r (e) mı r Be : ae 


Be Für das Fehlerverhältnis Al gilt mit (22) | war * BR t j 
2 ER | (ar al: IE Emsumtr) Rz a I 

dx SR k - 7 N FR r 2% . i 

ER | PD) I Im B77) 5 Kr > SI 


Er er 
Nun ist entsprechend zu (21) bei eben Werten ff vyundv S ar: 
Br Zi Mina a = Tr z uk» 
Re EN also | Kr iR ge | 
20 RER Pal Fir n / fh 
m... 5 3 Sau, ar 
FR Ir u v=—n v’=—n PER 
RR RL Setzt man hier ?’ =» +0, so kommt ae 4 
2 | Ne) nn +2n rn 
Ber: MN Dunn > Kt Zei i,.) 
Be,: a v-—n 0=—2n a a 
Be: oder wegen (24,) ; 
BR JE Kt r, DE 1er.) rer PAGE 
Ye 0=—2n g=1.- ı ; 
| Das Fehlerverhältnis iu hängt unter anderem von der bereits vorhandenen Korrelation r, 
* der Ausgangsfehler x und von der ‚‚durch das Ausgleichen erzeugten Korrelation“ r, ab. 
Ein paar Sonderfälle sind zur Veranschaulichung nützlich. Für, =0?e= +1; +2;...) 
! geht (29,) in die ‚normale‘ Fehlerformel 
‚ 2 “ 
| &) EN BR. 
über. Bei vollständiger Korrelation r, =1 wird (z/u)? = A,&r, oder wegen (25) x = u. Esist 
also kein Ausgleichen möglich. Bei erheblicher negativer Korrelation unmittelbar benachbarter 
Ausgangsfehler u, und a4, wird abwechselnd r, < 0, r3>0, 7, <0, r,>0,... und damit j 
“ 4 
&) = A f rd 2 (Irze-ılFae-ı — rae | . | 
Wenn nun die Korrelation r, zwischen 1; Ya fx+. mit wachsendem Abstand abnimmt (was . 
der Regelfall sein wird), so sind die einzelnen Glieder der Summe positiv und damit (z/4)? < },. J 
Der Ausgleich ist also besser als im ‚„‚Normalfalle“ (29,). 4 
” Das zwischen den Fehlern ji, der ER Werte &; bestehende X-Maß wird 1 
| - fix Hgg N 
en Ei 
To = Il etw. ee ee a ee (30); | 


wobei der angehängte * die Verwechslung mit der im Abschnitt 7-erklärten X-Funktion r, aus- 1 
schließen soll. In der letzten Gleichung ersetzt man f%, und +. durch (22) und beachtet, 2 | 
daß jetzt gemäß (28) im Gegensatz zu (23,, ,) 4 
! 2 
I Mir Mktotv! = Tot - Ik 


ist. Nach einer kurzen Rechnung, die ähnliehswie vorhin verläuft und deshalb nicht im einzelnen 
wiederholt werden soll, erscheint #* schließlich in der Gestalt 


--2n 2n 
... n 
2 Tat tr Hr 2 (ur Free 
-x _ 9=—2n 5 y=l 
a TEE (31). 


2 ryr, 1.2 3 
-- vi» 


Be 


PER: 4 we BR Ä KR j! \ x / 
Ba.0 eh" Stange, Mehrfaches Ausgleichen einer fehlerhaften Punktreihe 2.128 


T 


9. Mehrfaches Glätten bei unabhängigen Ausgangsfehlern. Wendet man (31) auf 
das mehrfache Ausgleichen ursprünglich korrelationsfreier Meßfehler an, so wird das K-Maß 


r0+» zwischen den Fehlern der (0 +1)-mal geglätteten Werte x,» 
+2n N 
(a) 
| HEN 
„tl _ v=—an srkfe | 32 
N DON N 
. PER 
v=—2n 
it (+1) _ (+1) _ n 
mi ro = Yundarn, =0 für Je|>2n(c-+1). 


“ Bezeichnet man den Nenner in (32) kurz mit Ns, so erhält man als Summe $,;, aller 
r0+D den Ausdruck | | 
e 


8 a (041) 1 zn EC: () e 
Lo a We = No = N. 2, Ty = Totp: i $ 
0= 0 co v=—2Nn 0=— © 
Nun ist bei festem » 
+ © (0) + © (0) ; 
D3 Toty = ER I So; 
e=—° A=— 80 
‚also 
I TE Ic 
Ss ee = S 
> o+1 Ns Abrmat A v Ns 2 
oder us N, = 1 
- . o-+1 
BD ee 
N, N ar N 


eine Beziehung, die man zu Prüfungszwecken heranziehen wird. 

Die K-Maße lassen sich mit der Rekursionsformel (32) von Schritt zu Schritt berechnen. 
Für lineares und quadratisches Glätten mit je N = 5 Punkten (bzw. n = 2) ist das in der folgen- 
den Zahlentafel geschehen !). Jedes K-Maß Ton ') der Spalte (o +1) ergibt sich als normiertes 
„gewogenes Mittel‘ der2%» +1 benachbarten Werte 9, u ne, 2, nn 2„ der vorausgehenden 
Spalte o. Als „Gewichte“ dienen die r, und normiert wird so, daß r(” unabhängig von o den 
Wert 1 erhält. Mit wachsender Schrittzahl o verbreitert sich die K-Funktion re mehr und 
mehr, wobei die der Mitte benachbarten für die Fehlerfortpflanzung maßgebenden K-Maße sich 
dem Grenzwert 1 nähern. Für ‚‚große‘‘ Werte von o wird also mit 25 


Neue I Beh we N 


Tafel der Korrelationsmaße ER = a: für N=5 


Lineares Ausgleichen Quadratisches Ausgleichen 
0 | 1 2 | 3 | 4 1 | 2 | 3 4 
e 
0 1 1 1 1 1 1 1 8 1 
(1 0 0,8 0,941 0,963 0,971 + 0,565 | + 0,709 | + 0,764 | -+ 0,794 
2 0,6 0,800 0,860 0,891 + 071| + 194 |-+ 287 ++ 355 
3 0,4 0,612 0,712 0,770 — 121 | — 086 | — 038 | + 007 
4 0,2 0,412 0,543 0,627 + 015 | — 061 | — 090 | — 095 
5 0 0,235 0,380 0,480 ) — 000 | — 0225| — 051 
6 0,118 0,243 0,344 + 010 | + 010 + 001 
ri 0,047 0,141 ‚0,231 — 002 | + 005 + 009 
8 0,012 0,072 0,143 — | — 001] + 002 
9 DR. 0,032 0,082 0 — 001 | — 00 
10 0,012 0,043 — | — 00 
11 0,003 0,021 — — 
12 0,001 0,009 — _ 
13°. d 0,003 0 a 
14 0,001 — 
No 1 3,400 4,120 4,360 4,490 1,678 1,848 1,910 1,942 
So 1 5 7,353 8,924 10,234 2,06 2,53 2,82 3,04 


11) Die gegebenen Zahlenwerte besitzen ‚‚nur Bechenschiebergenauigkeit‘“, 


Stange, Mohn 


a Das Verhältnis zweier a nei 
BERG Km 1 . S 2 2 H ie get n - 3 DR 
d.h. das Glätten wird mit wachsender Schrittzahl wirkungslos, was bereits anscha 
0 Bild 3 hervorgeht. | IR | Er 
10. Ausgleichen einer fortlaufend gegebenen Kurve bei korrelativer } 
> 0... Ausgangsfehler. In Bild 4 soll der geglättete Wert =, in der Mitte des festen Bei 
eh W = — B<ß<+Bau N=2n+ 


Re punkten 2_,... 2... Zn 

5 Pr werden, wobei wir uns vorbehalten, 
2 . die Punktzahl n StR ae ar 
Ban: schränkt wachsen lassen. Wir g- 
en ! langen dann zum Ausgleichen einer rt 
B; nicht mehr durch Einzelpunkte son- 


dern fortlaufend gegebenen Kurve), 
Daß dannbenachbarte Meßfehlernicht 

mehr unabhängig voneinander sen 

können, ist wegen der bei jedem Meß- 

gerät vorhandenen Trägheit einleuch- 

tend. Die Größe von B ergibt sich 

OR aus der Forderung, daß der mathema- 

Bild4. Zum Ausgleichen einer fortlaufend gegebenen Kurve. tische Fehler oder die Fälschung des 

| geglätteten Wertes innerhalb erträg- 

licher Grenzen bleiben muß. Die zwischen benachbarten Fehlern bestehende Korrelation kenn- 

zeichnen wir durch die gerade Funktion r (ß), welche der Forderung genügt 


| r(ß) #0 für |Bl<D, | 
mit Ausnahme endlich vieler Nullstellen, fer ren RE N 
"(d)=0. für |Bl>b 


Die Größe 5b nennen wir Korrelationsbreite, und 


+b b 
Fr | Bar=z (rmar Er eier 3 (35,) 
—b 0 


heißt. das mittlere X-Maß. rund d sind zwei das Meßverfahren und das Meßgerät kennzeich- 
nende Größen. Da 5 normalerweise ‚klein‘ gegen Bist, so dürfen wir voraussetzen, daß der zu 
glättende Bereich 2B ein geradzahliges Vielfaches der X-Breite 5 ist, 


B=n,d, N. = ganze. Zahl, nur un (36). 
Ist nun zunächst beim Ausgleichen n— n, < n,, so ist in (29,) die Summe für alle oe #0 
leer und es gilt die ‚‚normale“ Fehlerformel (29,). Überschreitet jedoch » die „‚kritische‘‘ Punkt- 


zahl n,, so gilt das Fehlergesetz (29,), das wir im Hinblick auf den Grenzübergang n— © 
zunächst umgestalten. Wir schreiben 


ak vB 


MH 2B AV us. "e n j 
und setzen 1 
SER B B | | 
| ß — AR bzw. AB=—4e= u EEE ne Dae- (37). | 
Dann wird ; 4 
A 2R A, te e \ 
= DAB. 
Lassen wir jetzt n über alle Grenzen wachsen, so strebt 2nı, gegen den Grenzwert ®) 
i Be I u u ED 
(2n A,» E DIET — De 1. 2 | 


Dabei ist R= 2p bzw. KR’ =2p--1 die Ordnung der beim Ausgleichen benutzten Ersatz- 
funktion. Für die hier in Betracht kommenden einfachsten Fälle des linearen (p=0) und 


22, Willers,'a.a. 0,78. 254, 
13) Vgl. Fußnote ®). 


eh nun erh 11 Aal na 
\ 
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quadratischen Ausgleichs (p = 1) gilt insbesondere 


9 
@rd)gz=1. und, Ar); > Be AR ee ee (38). 
Das Fehlergesetz im Grenzfalle wird also | 
& . B2+42B | 
.(ü\® in d)o e 
Mens wa 


p=—-2B 


oder, da sowohl r (ß) als auch r (ß) gerade Funktionen von ß sind und r (8) für || > d identisch | 
verschwindet, 


| | (BAR. = Di ae en. 


Von $=0 bis = 2B fällt r (ß) vom Werte r (0) = 1 auf den Wert 7 (2B) =0 ab. Im Inte- 
grationsbereich O<Pß<b ist BeseH b<2B nahezu 


r(p) 1 Re NR an (39,). 
Das Fehlerverhältnis wird also mit (35,) näherungsweise 


Er (P)dP = (2nA,) NERZTSE, a (395): 
P ; 


Während der Fehler z der ausgeglichenen Werte bei korrelativer Unabhängigkeit der 
Ausgangsfehler «u mit wachsender Punktzahl n beliebig klein wird, strebt er bei korrelativer 
Verknüpfung der Ausgangsfehler einem im allgemeinen von Null verschiedenen Grenzwert u, 


- zu, dessen Quadrat u, unter anderem dem mittleren X-Maß ry, und dem Quotienten der Bereiche 


b/B verhältnisgleich ist. 
Für den bei der „kritischen“ Punktzahl n, erreichten Fehler 4, gilt nach (29,) 


u z et (an Ao)x ae b 
(2), — ans: — (2n Ao)s SB a Ta BE A ET a ON (40). 
Führt man 4, in (39,) ein, so kommt 

R (ANA) as 

lo — (en ar J Arm Ux ee re (41,). 


“Die erste Wurzel hat beim linearen Ausgleichen den Wert 


an Fl AN al 
Y). en a ne 


und beim quadratischen Glätten- gilt 


DR am DO FHOM+N_,, 1 er! ” 
a Sa (au an, Ma N go mE 


B j 2 i A 
Dan, = y® normalerweise eine ‚große‘ Zahl ist, so darf man die Wurzel in beiden Fällen mit 


guter Annäherung durch den Wert 1 ersetzen. 


RE TE NE (415). 


Dann kommt 


Ist das mittlere X-Maß ry > x so wird das Ergebnis des Ausgleichens, verschlech- 
tert, wenn man das Verfahren über den ‚‚kritischen“ Wert n=n, hinaus fortsetzt. 
Für ru= 5 wird der günstige Einfluß der wachsenden Punktzahl durch die sich mehr 
und mehr ungünstig auswirkende Korrelation gerade ausgeglichen. Das Verfahren läuft für 
n >.n, „leer”. Fürry < “ ist zwar hoch eine Verbesserung der Streuung möglich. Jedoch 


wird sich der erforderliche Mehraufwand an Rechenarbeit im allgemeinen nicht lohnen. Eine 


ers 
‚Ein Beispiel mag das Vorausgehende noch 
man gemäß co) und (37) im Grenzfalle die 


EIER Ss lim 7, — lim ( Ber 


oder 


Ta ER N R: 
ale 128 6 B 93B Er re \ R 
| RE: og RR 

Für den bei der kritischen Punktzahl n, = „ erreichten Fehler ., gilt x: | 
2 Te I, \ a: = RI == Ze y Br ( u iR & > % a VE 3 (42). h "* N | 
Bi BR, er re “zB : 2B BEN Fr Er 
B: iR, » B 1 b PER 
E Der geringe hier zwischen i und x, auftretende relative Unterschied = ‚ der tatsächlich 
A: | 


ohne Bedeutung ist, wird in der Näherungsformel (41,) durch die Vernachlässigungen (39,) und 


b : 
H (41,) verwischt. Für das Bereichverhältnis 3” 0,2 ist der Verlauf von # über der Punkt- f 
Be zahl n in Bild 5 dargestellt. Bis zum „kritischen“ Punkt P,(n,) sinkt der Fehler # „normal“. 


mein 


a 


= N Er Bert 


7 5 h 70 DB i 
Bild 5. Beim linearen Ausgleichen mit wachsender Punktzahl n in einem festen Bereich B sinkt der Fehler a der geglätteten 
Werte bei „dachförmiger‘ Korrelation der Ausgangsfehler « nur bis zum „kritischen“ Punkt n = n,. 


Von da ab schwingt der tatsächliche Verlauf mit geringer Amplitude um den Grenzwert (ul) » 
Die Schwingungen veranschaulichen gewissermaßen den Kampf zwischen dem günstigen Ein- f 
fluß der wachsenden Punktzahl und dem hemmenden Einfluß der Korrelation, die sich mit j 
wachsender Punktzahl auf mehr und mehr Punkte ausbreitet, . 
Ein sinnvolles Ausgleichen fortlaufend gegebener Kurven ist also nur dann möglich, wenn 
man neben der Größenordnung der Ausgangsfehler u auch die Größenordnung ihrer Korrelations- 4 
breite b kennt. Dabei ist es zweckmäßig, die in einem festen Bereich B benutzte Punktzahl n 


nicht größer als den „kritischen“ Wert n, -7 zu wählen. 


14) Welcher der drei eben genannten Fälle bei fortlaufend 
nicht untersucht werden, da das Beobachtungsmaterial, an dem 
durch Kriegseinwirkung verloren ging. 


Eingegangen: 25. Mai 1948. 


gegebenen langen Meßreihen eintritt, konnte 
die Untersuchung durchgeführt werden sollte, 


"Zur Maximalkorrelation. 


Unter dem Titel „Zur Maximalkorrelation‘‘ haben 
G.Friede und H.Münzner in dieser Zeit- 
schrift?) den folgenden Satz bewiesen: 

Hinreichend dafür, daß die Maximalkorrelation K 
den Wert 1 annimmt, ist, daßin der Korrelationsmatrix 
alle Elemente, in denen sich n beliebige Spalten mit m 
beliebigen Zeilen schneiden, und alle Elemente, in 
denen sich die restlichen k—r Spalten mit den rest- 
lichen I—m Zeilen schneiden, gleich Null sind. 

Wir werden im folgenden zeigen, daß diese Be- 
dingung auch notwendig ist, und daß ihre Notwendig- 
keit und ihre Hinreichen in ganz einfacher Weise direkt 
aus der Definition von K folgen, ohne daß die Ver- 
wendung der Determinante F(z) erforderlich ist, die 
in der zitierten Arbeit herangezogen wurde. 

Da der Wert von K definitionsgemäß invariant ist 
gegenüber Permutationen der Zeilen und Spalten der 
Korrelationsmatrix, läßt sich die Friede-Münz- 
nersche Bedingung auch so ausdrücken, daß die 
Korrelationsmatrix N durch Umordnen der Zeilen und 
Spalten in die folgende Gestalt gebracht werden kann: 


(ea) 


wo in den Rechtecken R, und R, beliebige Zahlen 
stehen. 

Eine solche Matrix wollen wirreduziertnennen. 
N heißt reduzibel, wenn sie reduziert werden 
kann. 

In Erweiterung des Satzes von Friede und 
Münzner:gilt nun: 

Notwendig und hinreichend dafür, 
daß K=l ist, ist die Reduzibilität 
vonN. 

In dieser Gestalt ist der Satz ganz analog dem Satz 
von H. Gebelein?) über kontinuierliche Ver- 
teilung, wo übrigens unter Vermeidung des in seinem 
Lehrbuch $S.381 stehenden Versehens ausdrücklich 
bemerkt wird, daß bei K=1 für den stochastischen 
Zusammenhang ‚noch ein gewisser Spielraum- ent- 
steht, den K nicht zu erkennen gibt‘, 

Zunächst ist nun das von FriedeundMünzner 
bewiesene Hinreichen vollkommen evident. Setzen 
wir nämlich die Merkmalswerte für ‚alle Zeilen und 
Spalten, die durch R, laufen, einerseits und für alle 


1) Z. angew. Math. Mech. Bd. 28, 1948, $.158—160. Der Ein- 
fachheit halber übernehmen wir die Bezeichnungen dieser Arbeitohne 


neue Erklärung. 
2) Z. angew. Math. Mech. Bd.21 (1941), 8.371. 


i 


KLEINE MITTEILUNGEN 


durch R, laufenden Zeilen. und Spalten andererseits je \ 
untereinander gleich, so erhalten wir eine neue Korre- 
ationsmatrix N’, die die Gestalt 
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mit dem Korrelationskoeffizienten K’ = 1hat. Wegen 
1> K> K’ ist also auch K =1. 
Ist umgekehrt X =1, so folgt: faus 


ID ea N; 1 


vi 5 ' 
wegen der gleichzeitig erfüllten Normierungsbedin- 
gungen ' 


Dre N4=- N eWNg=N 
%) 65 
sofort: f(zi) = g(yj) für alle N,; 0. 

Ohne Einschränkung der Allgemeinheit können wir 
annehmen, daß N,,0 ist. Wir ordnen nun die 
Zeilen und Spalten ‚von N so um, daß diejenigen 
Zeilen und Spalten als die ersten erscheinen, bei denen 
fe) = gu) =) =g4n) ist. Wegen der Zen. 
trierungsbedingungen 


De) Ny = Da N;=0 
07 4,5 


bleiben dabei sicher noch Zeilen und: Spalten übrig 

In den Schnittpunkten dieser restlichen Zeilen mit 
den’ersten Spalten ist f(x) g(yj), also N; =0. 
Entsprechendes gilt für den Schnitt der restlichen 
Spalten mit den ersten Zeilen. Damit ist also gezeigt, 
daß N eine reduzierte Gestalt angenommen hat. 

Aus dem Dargelegten. folgt insbesondere, daß für 
irreduzible N stets K?<(1 gilt. Hieran kann man 
noch eine einfache Bemerkung knüpfen. Da für die 
reduzierte Gestalt von N offenbar 


F@A)=P,@)-F,@) 


gilt, wo F,(z) und F, (2) resp. zu den Rechtecken R, 
und R, gehören, ist die Anzahl der Lösungen z=]1 
der Gleichung F (2) = 0 gleich der Anzahl der Recht- 
ecke, in die N zerfällt, wenn man so lange weiter aus- 
reduziert, bis nur noch irreduzible Rechtecke übrig 
bleiben. Hieraus folgt, daß diese Anzahl der irre- 
duziblen Bestandteile unabhängig ist von der Art 
ihrer Gewinnung, eine von der Algebra her geläufige 
Tatsache. 


Haltingen (Kreis Lörrach). H.Richter. 


BUCHBESPRECHUNGEN 


M.H. Sassenfeld, H. F. A. Tschunko: Mathema- 
tische Tafeln für Mathemätiker, 
Naturwissenschaftler, Ingenieure. 
VIII+ 34 S. Walldorf bei Heidelberg 1949, Fr. 
Lamade-Verlag. 5,75. DM. 

Der vorliegende erste Teil mathematischer Tafeln 
enthält zunächst die üblichen elementaren Funk- 
tionen (Potenzen, Wurzeln, Kehrwerte, Kreis-, 
Hyperbel-, Exponentialfunktionen, Briggsche und 
natürliche Logarithmen sowie auch zyklometrische 
Funktionen und Antilogarithmen zwecks Vermeidung 
inverser Interpolation), und zwar vorwiegend fünf- 
stellig bis auf die vierstelligen Briggschen Logarith- 
men nebst Winkelfunktionen nach Grad und Neugrad. 
Weiter werden Umrechnungstafeln für Winkeltei- 
lungen, Tafeln für Kreis- und Kugeldaten, eine Tafel 
der Primfaktoren und der sechsstelligen Logarith- 
men der Primzahlen fürn = 1—2000 sowie Vielfache, 


Potenzen, Logarithmen von z, M, e, oe, Potenzen und 
Fakultäten von n, Binomialkoeffizienten u. a. (zehn- 
bis zwanzigstellig) gebracht. — Bei den Werten tgh x 
hätte man vielleicht für größeres Argument (z > 2) 
durch Übergang auf 1—tghz auf gleichem Raum 
eine noch größere Stellenzahl wiedergeben können. 
Eine Neuerung besteht darin, daß durch geschickte 
Satzanordnung in Verbindung mit DIN A 4-Format 
101 Argumentzeilen auf einer Tafelseite untergebracht 
werden, wodurch auf einer Seite vier gewöhnliche 
Tafelseiten Platz finden, ohne daß die Übersiehtlich- 
keit leidet. Der Argumentschritt ist so gewählt, daß 
fast durchweg linear interpoliert werden kann. Was 
die Tafeln aber für den wissenschaftlichen Rechner 
besonders wertvoll macht, ist die Aufnahme einiger 
meist zehnstelliger Zusatztafeln, mit deren Hilfe man 
sich bei mäßigem Arbeitsaufwand (zum Teil mit dem 
Rechenschieber) zu fast allen aufgeführten Funk- 


128 


tionen auch vielstellige (meist zehn-, oft aueh noch“ 


mehrstellige) Werte verschaffen kann, ohne auf schwer 
zugängliche umfangreiche Tafelwerke angewiesen Zu 
seln. 
Anleitung, nebst Beispielen und Fehlerabschätzungen 
gibt hierzu'auch dem mathematisch weniger Geübten 
alle erforderlichen Erläuterungen. j 

“ Man möchte dem sorgfältig angelegten und gut 
ausgestatteten Tafelband die verdiente Verbreitung 
wünschen und darf auf die angekündigten Fort- 
setzungen (logarithmische und logarithmisch-trigono- 
metrische Funktionen, natürliche trigonometrische 
Funktionen, spezielle Funktionen) gespannt sein. 

Darmstadt. R. Zurmühl. 


Dr.-Ing. Fr. Boönjakoviöd (o. Prof. an der Universität 
Zagreb). Technische Thermodynamik (Wärmelehre und 
Wärmewirtschaft in Einzeldarstellung. Herausgegeben 
von Dr.-Ing. W. Pauer, Band 11, 1. Teil). Dritte un- 
veränderte Aufl. XVI+ 327 8. mit 249 Abbildungen und 
3 Tafeln. Dresden und Leipzig 1948. Verlag Theodor 
Steinkopff. Preis: geb. 18,— DM; brosch. 17,— DM. 

Die erste, 1935 erschienene Auflage dieses allgemein 
als vorzüglich anerkannten Werkes wurde in dieser 
Zeitschrift, Bd.16 (1936), 8.126 von Trefftz be- 
sprochen. Eine 2. Auflage erschien 1943. In ihr sind 
die Abschnitte über Verbrennung und Vergasung aus 
didaktischen Gründen aus dem 2. in diesen 1. Band 
übernommen, und es ist ein ergänzender Abschnitt über 
Wärmeübertragung hinzugefügt worden. Die vor- 
liegende 3. Auflage wurde aus zeitbedingten Gründen 
als anastatischer Neudruck der 2. Auflage hergestellt. 

Was Trefftz lobend von der 1. Auflage sagt, daß 
die Darstellungsweise sich durch Anschaulichkeit und 
strenge Deduktion auszeichne, das gilt auch von den 
neu hinzugefügten Abschnitten. Man muß dem Verlag 
dankbar sein, daß er das Buch durch diesen Neudruck 
den deutschen Wissenschaftlern und Ingenieuren wieder 
zugängig gemacht hat. 

Dresden. Willers. 

Dr. Heinrich Blasius (Oberstudienrat a. d. Ing.- 
Schule Hamburg),Mechanik, physikalische Grund- 
lagen vom technischen Standpunkt, erster Teil Sta- 
tik ‚4. Aufl., VIII + 247 S. mit 260 Fig. u. 113 Aufg. 
Hamburg 1948. Verlag Eckardt u. Messtorff, brosch, 
DM 9,—. 

Der erste Band ‚‚Statik‘ der Buchreihe des Verfassers 
über ‚‚Mechanik, physikalische Grundlagen vom tech- 
nischen Standpunkt‘ hat nunmehr seine vierte Auflage 
erlebt, Infolge der reichhaltigen Auswahl der aufge- 
führten Beispiele wird der Leser zur selbständigen Lö- 
sung der Aufgaben angeregt und so zwanglos mit den 
einfachsten Grundprinzipien des Fachgebietes vertraut 
gemacht. So besitzt das Buch einen verhältnismäßig 
hohen pädagogischen Wirkungsgrad, wenn auch die 
Art der Vermittlung der Grundzusammenhänge etwas 
abwegig ist und zum weitergehenden Studium der Me- 


ine den Tafeln vorangestellte übersichtliche _ 


inzi weh eniger wichtige 
Tafeln der er des 
ferner Tafeln der Laguerr 


die Vandrey in dieser Zeitschrift 
und die ebenfalls in der Z. angew. Math. 
schienenen Tafeln der unvollständigen r | 
und Weberschen Funktionen von P. undE.Brauer. 
Die Dinnikschen Tafeln für Jn);(iy) und J_nfs in sind 
durch neue von Kerridge berechnete für Jn/,(4y) 
und HU), (iy) (m=1; 2) ersetzt. Weiter sind Formeln _ 
und Figuren zung Gebrauch der 
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Debyeschen, gi = I 
konvergenten Reihen der Zylinderfunktionen RUE 
komplexen Argument und Index ei die Tafeln Y 
für das Fehlerintegral erweitert und ein] weisvon 
sieben- und mehrstelligen Logarithmentafeli an 


gehängt. i ET 
Fortgelassen sind dagegen unter anderem & 
über der dritten Auflage die Abschnitte vi. 14 
(Spezielle Weierstraßsche ktionen für g, = 0 und 

„= ]) und VIIL,8. (Die Sommerfeldsche und die 
Funktionen.) ae 


oissonsche Integraldarstellungen der 

Die vor allem für die numerische Auswertung bei 
Anwendungen der Mathematik auf physikalisch und 
technische Probleme unentbehrlichen Tafelwerke wer- 
den nach wie vor von Mathematikern, Physikern und 
wissenschaftlich arbeitenden Ingenieuren viel benutzt 
werden. Alle Benutzer werden Herrn Emde für 
die Mühe dankbar sein, die er auf die immer bessere 
Ausgestaltung seines Werkes verwendet. 


Dresden. Willers. 


Die besprochenen und angezeigten Bücher sind durch den Buchhandel zu beziehen \ 


NACHRICHTEN 


Leipzig: Mit Wirkung vom 1. Oktober wurde der 
frühere o. Prof. an der Universität Königsberg, Dr. 
ErichKähler (zuletztin Hamburg) zum Professor 
mit Lehrstuhl für Mathematik in der philosophischen 
Fakultät der Universität Leipzig ernannt. 


München: Prof. Dr. R. Sauer (früher Aachen) 
wurde zum 1. Oktober 1948 als 0. Professor für höhere 
Mathematik und analytische Mechanik auf den Lehr- 


Verantwortlich für den Inhalt: 


stuhl von Prof. Baldus, vorher v. Dyck an die 
Terhnische Hochschule München berufen. 


‚ Dozent Dr. W.Damköhler hat einen Rufan 
die Universität Tucumän (Argentinien) angenommen. 


Frankfurta.M.: Dr. Bernhard Mrowka 
wurde mit Wirkung zum 29. 6. 1948 die venia legendi 
für theoretische Physik verliehen. 


Prof. Dr. Fr. A. Willers, Dresden; für den Verlag: H. Kaesser, Berlin, 


Verlag: Akademie-Verlag GmbH., Berlin, NW 7, Schiffbauerdamm 19, Fernsprecher: 42 6412 und 426918. Postscheckkonto: 


Berlin 35 021. 
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